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Abstract. We investigate the convergence of a nonlinear approximation method introduced by 
Ammar et ah (J. Non-Newtonian Fluid Mech. 139:153-176, 2006) for the numerical solution of 
high-dimensional Fokker— Planck equations featuring in Navier-Stokes-Fokker— Planck systems 
that arise in kinetic models of dilute polymers. In the case of Poisson's equation on a rectan- 
gular domain in R-^, subject to a homogeneous Dirichlet boundary condition, the mathematical 
analysis of the algorithm was carried out recently by Le Bris, Lelievre and Maday (Const. 
Approx. 30:621-651, 2009), by exploiting its connection to greedy algorithms from nonlinear 
approximation theory, explored, for example, by DeVore and Temlyakov (Adv. Comput. Math. 
5:173—187, 1996); hence, the variational version of the algorithm, based on the minimization of a 
sequence of Dirichlet energies, was shown to converge. Here, we extend the convergence analysis 
of the pure greedy and orthogonal greedy algorithms considered by Le Bris et al. to a technically 
more complicated situation, where the Laplace operator is replaced by an Ornstein— Uhlenbeck 
operator of the kind that appears in Fokker— Planck equations that arise in bead-spring chain 
type kinetic polymer models with finitely extensible nonlinear elastic potentials, posed on a 
high-dimensional Cartesian product configuration space D = Di X ■ ■ • X Djv contained in ] 
where each set Di, i = 1, N , is a, bounded open ball in R"*, d = 2,3. 



1. Introduction 

High-dimensional partial differential equations are ubiquitous in mathematical models in sci- 
ence, engineering and finance. They arise in a number of areas, including, for example, kinetic 
theory, molecular dynamics, quantum mechanics, and uncertainty quantification based on poly- 
nomial chaos expansions, to name only a few. 

The purpose of the present paper is to explore the convergence of a numerical algorithm that 
was recently proposed in the engineering literature in a succession of papers by Ammar, Mokdad, 
Chinesta, Keunings and collaborators |AMOK06[ IAMOKOTI lAND+lOl ICjACClOl ICALKllj . for 
the numerical solution of high-dimensional Fokker-Planck equations in kinetic models of polymeric 
fluids under the names Separated Representation and Proper Generalized Decomposition. A variant 
with a discretization based on spectral methods instead of the finite element methods preferred by 
Ammar et al. was presented by Leonenko and Phillips |LP09) . A similar method was considered 
independently by Nouy |Nou07l INou08| and Nouy & Le Maitre [NLMOQj under the name Power 
type Generalized Spectral Decomposition, for the numerical solution of stochastic partial differential 
equations, although the historical roots of the technique can be traced back to the work of Schmidt 
[Sch07| . Ammar et al. and Nouy report that the algorithm performs well in numerical experiments 
and comment that it extends to a large variety of partial differential equations. 

In the simplified mathematical setting of Poisson's equation —Au = f posed on the rectan- 
gular domain n = x fly, where and Vly are bounded open subintervals of R, subject to a 
homogeneous Dirichlet boundary condition on dfl, the convergence of the algorithm was shown 
in a recent paper by Le Bris, Lelievre and Maday |LBLM09) . by drawing on connections with 
greedy algorithms from nonlinear approximation theory (cf. DeVore and Temlyakov |DT96) ). In 
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[LBLM09] . the solution was represented as a sum 



u{x,y) = ^r„(a:) Sn{y) 



(1.1) 



n>l 



by iteratively determining functions x € '-^ fnix) and y £ f2y i— )■ s„(y), n > 1, such that for 
all n, the product (x,?/) € i— >■ r„(a;) s„(y) is the best approximation in the norm of the Sobolev 
space Hj(il) to the solution {x,y) e H> v{x,y) of the Poisson equation 



subject to a homogeneous Dirichlet boundary condition, in terms of a single function of the 
factorized form r{x) s{y); Le Bris et al. thus show that it is possible to give a sound mathematical 
basis to the algorithm proposed by Ammar et al., provided that one considers a variational form 
of the approach that manipulates minimizers of Dirichlet energies instead of stationary points to 
the associated Euler-Lagrange equations (in the follow-up paper |CELllj by Cances, Ehrlacher 
and Lelievre it was further shown that one can also work with local — yet still energy-decreasing — 
minimizers provided that one stays within the two- fold tensor product setting of (jl.ip ). In order 
to reformulate the approach in such a variational setting, the arguments in |LBLM09] crucially 
rely on the fact that the Laplace operator is self-adjoint, and as noted by the authors of [LBLM09] . 
the analysis does not apply exactly to the actual implementation of the method as described in the 
papers by Ammar et al., where stationary points of the Euler-Lagrange equations associated with 
the Dirichlet energies are computed instead. Indeed, since minimizers of Dirichlet energies in the 
approach of Le Bris et al. on the one hand and stationary points of the associated Euler-Lagrange 
equations in the approach of Ammar et al. on the other are each sought in nonlinear manifolds 
embedded in a Sobolev space, rather than over the entire Sobolev space (which is a normed 
linear space), the two approaches are not equivalent. The authors of [LBLM09J also comment 
that: "Likewise, it is unclear to us how to provide a mathematical foundation of the approach 
for nonvariational situations, such as an equation involving a differential operator that is not self- 
adjoint." This latter remark is particularly pertinent in the context of Fokker-Planck equations 
for kinetic bead-spring chain models for dilute polymers, of the kind considered by Ammar et al., 
where the differential operator in configuration space featuring in the Fokker-Planck equation, a 
generalized Ornstein-Uhlenbeck operator, is a non-self-adjoint elliptic operator with a drift term 
that involves an unbounded potential. 

It is this last point that the present paper is aimed at addressing: we shall be concerned with the 
numerical approximation of high-dimensional Fokker-Planck equations that arise in bead-spring 
chain type kinetic models of dilute polymers on the Cartesian product domain x D, where 
n C M'', d = 2,3, is the physical (flow) domain, and the configuration space D is the A^-fold 
Cartesian product Di of sets Di C K"^, ? = 1, . . . , iV, > 2, each of which is a bounded open 
ball in W^. Here, denotes the numbe r of springs connecting, in a linear fashion, the A^-|-l beads 
in the bead-spring chain model (cf. Fig.O). Proceeding as in [BSOTllBSOSllBSMllBSllallBSllbj . 
we rewrite the Ornstein-Uhlenbeck operator, a non-self-adjoint elliptic operator with respect to the 
configuration space variable q featuring in the Fokker-Planck equation whose drift term contains 
an unbounded potential, as a degenerate, but now self-adjoint, elliptic operator on a Maxwellian- 
weighted Sobolev space. We then perform a nonlinear approximation of the analytical solution 
ijj: {qi,. . . , <7jy) e D V(9ii ■ • ■ > Qn) to this high-dimensional degenerate elliptic boundary-value 
problem on the appropriate Maxwellian-weighted Sobolev space by separated representations of 
the form 



where the factors V'fc j ^ = 1; ■ • • i are defined on the d-dimensional domain Di, i — 1, . . . ,N. 
Instead of being selected from an a priori fixed set, the factors ipf^ , i — 1, . . . ,N, are obtained, 
A^ at a time, for each G {1, . . . , K}, as the best approximation (in a sense to be made precise in 




K N 




k=l i=l 
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Section 0) among all possible such factors. The (potentially large) number of terms K is likewise 
not fixed in advance, but depends on a termination criterion. 

The paper is structured as follows. After introducing our notational conventions and formu- 
lating briefly an alternating direction scheme that separates, by a fractional step method, the 
full Fokker-Planck equation into a low-dimensional physical space part and a high-dimensional 
configuration space part, we will concentrate on the latter problem. The central difficulty in the 
numerical solution of the configuration space problem is the presence of the high-dimensional 
Ornstein-Uhlenbeck operator, a non-self-adjoint elliptic operator whose drift term contains an un- 
bounded potential. In Section we show that the configuration space problem can be restated, in 
a Maxwellian- weighted Sobolev space, as the weak formulation of a symmetric degenerate elliptic 
boundary- value problem on the high-dimensional configuration space D. Section |3| is devoted to 
the description of a separated representation strategy for the problem, in the spirit of Le Bris 
et al. |LBLM09] . Followin g_|LBLM09| . we consider a pure greedy algorithm and an orthogonal 
greedy algorithm. Section concentrates on the convergence of the two algorithms. We shall 
characterize the convergence rates of the two greedy algorithms by invoking abstract convergence 
results due to DeVore and Temlyakov |DT96| . In Section [HI, we give explicit necessary and suffi- 
cient conditions, in terms of Maxwellian-weighted Sobolev spaces, for membership of the space of 
DeVore and Temlyakov in the case of our degenerate elliptic problem. In Section 0, we provide 
some conclusions and possible directions for further work. 

At an abstract level, our convergence proof follows the arguments in |LBLM09] : however, the 
verification of certain key properties of the function spaces involved, on the one hand, and the 
characterization of verifiable sufficient conditions under which the predicted convergence rates of 
the two greedy algorithms considered are observed, on the other, for the high-dimensional degener- 
ate elliptic problem studied herein are considerably more complicated than in the case of Poisson's 
equation studied in |LBLM09] . The former is mostly based on tensorizing the corresponding re- 
sults for the function spaces associated with the single-spring case (i.e., the dumbbell) and the 
latter relies on shift-theorems for degenerate elliptic operators in Maxwellian-weighted Sobolev 
spaces and delicate results from the spectral theory of self-adjoint degenerate elliptic operators, 
which we were unable to find in the literature; these are described in Section and Appendix 0, 
respectively. Appendices O and |B] collect a number of technical results that are used throughout 
the paper. 

1.1. Notation. We denote by [k] the integer interval {i E N: I < i < k}. We shall denote 
sequences and arrangements of elements indexed by indices i in an index set I by {ai)^^.j-. 

We shaU write q = (q^, . . . , q^) € I?i x • • • x = X^gj^v] Di =: D. Given N real-valued 
functions /i, each defined on the corresponding set Di, we denote by (2)ig[Af] /* their tensor product] 
i.e., the function 

g e D ^ [] JMi)- 

We extend this notation in three ways. Firstly, as the tensor-product operation is order-dependent, 
we will use subscripts on the (g) and the signs to denote where on q G D the function, or func- 
tions, following them act; e.g., 0iG[Ar]\{j} ®j fj evaluated on q e D is fj{qj) nje[Af]\{j} fMi)- 
Secondly, we will use the same notation for the sets resulting from the tensor products of members 
of function spaces: suppose that Fi is a nonempty set of real- valued functions defined on Di, 
i £ [N]; we then write {S)ig[Ar] Fi :— {{^jgj^j fi'- fi € Fi, i e [N]}. Thirdly, if exactly one of the 
factors is vector- valued, the products involving it at the time of evaluation must be interpreted as 
scalar-vector products implying that the resulting tensor product will be vector-valued too. 

The symbol (s will stand for the compact embedding relation. The support of a real-valued 
function / will be denoted by supp(/). 

Given a measurable and almost everywhere positive real- valued function w defined on an open 
set C M"; i.e., a weight, we denote by L^(i?) the Lebesgue space of square-integrable functions 
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Figure 1.1. Bead-spring chain with TV springs and iV+1 beads. Adapted from 
Figure 11.4-1 of (BCAHSTj 



with respect to the weight w, equipped with its usual norm, 



\ip\'^w 



1/2 



We also define the w-weighted Sobolev space H™(i?) and its norm IHIh^cb) 

R[:{E) e Ll{E) n Ll,{E): d^^ E LI{E), \a\ < m} , 



|2 



1/2 



V^eH™(i?). 



^ |a|<m 

We shall suppose henceforth that 17 is a bounded open set in M."^ with a sufficiently regular (say, 
Lipschitz continuous) boundary, and denote by and rig. the unit outward normal vector defined 
(a.e. with respect to the surface measure) on dil and dDi, i e [N], respectively. 

1.2. Fokker— Planck equation. The spring forces in the model are given by functions Fi : Di — > 
R'^, which have the form Fi{p) = U',{\ \p\^)p, p £ D, :^ B{0,Vh) C K'^, b, > 0, i e [N], and 
the Ui: [0, 6,;/2) — ^ M, the spring potentials, are such that Ui{s) — > +oo as s — > fei/2_. It follows 
that Fi{p) = —Fi{—p) for all p E Di. Typical examples include the FENE (Finitely Extensible 
Nonlinear Elastic) model |War72] with 

C/,(s) = -|lnfl-^') and = — (1.2) 



where 6^ > is a parameter, and Cohen's Fade approximant to the Inverse Langevin (CFAIL) 
model [CohQlj with 



and 



l-k:lV(360 
^ iKQi) - — ; — — 9i 



(1.3) 



1 - W^V/h 

where 6^ > is again a parameter. We note in passing that both of these force laws are approxi- 
mations to the Inverse Langevin force law [KG42j 



where the Langevin function L is defined by L{t) :~ coth(i) — 1/t on [0,oo). As L is strictly 
monotonic increasing on [0,oo) and tends to 1 as its argument tends to oo, it follows that the 
function IgJ G [OjV^) ^~^(|9il/v^) € [O7O0) is strictly monotonic increasing, with a vertical 
asymptote at IgJ = ^/bi. 

Remark 1. An important spring force model, which is excluded from our considerations, is the 
simple Hookean model described by 



D, 



U^{s) 



and 



However, in many practically relevant flow regimes the physically unrealistic allowance of the 
Hookean model for indefinitely extended springs outweighs its mathematical convenience. 
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The Fokker-Planck equation under consideration for the probabihty density function -0 has the 
following form (cf. [B507l IBSOSl IB509l IBSllbllBmi] ): 

dib ^ r 1 ^ 

^ + div,(MV) + divq, (V,M)q,i^ - + Vq^. V) 



4Wi 



4Wi(A^°!-l) ^"^'' (a:,g,t)eOxDx(0,r], (1.4a) 



with initial and no-flux boundary conditions 

V'(-,-,0) = Vo, 



4Wi(A^+l) 



VccV ■ "a, = 0, 



and 



1 ^ 



(a;,g)er2xD, (1.4b) 
(a;,q,i) e 5f) X D X (0,r], (1.4c) 

•nq^=0, i e [A^], (a;,q,t) e r2xaDx(0,T]. (1.4d) 



Here, u : x [0,r] ^ R** is the flow velocity, Wi := XUo/Lo is the (nondimensional) Weissenberg 
number, Iq is the characteristic length-scale of a spring, A is the characteristic relaxation time of 
a spring and Lq and Uq are the characteristic macroscopic length and velocity, respectively (thus, 
Wi is the ratio of the microscopic to macroscopic time scales). The matrix A — (^,j )^ ^.^jj^j is 
symmetric and positive definite; we denote the smallest eigenvalue of A by Amin- 

We remark that the boundary condition (jl.4dp is an ensemble of N boundary conditions, which 
collectively account for the full {Nd— l)-dimensional measure of dD. 

We define the partial Maxwellians Mi and the (full) Maxwellian M by 



M,(p) :=Zriexp(-C/,(i|p|2)), p e D,, i e [N]; 



N 



Miq):=Y[M,{q,), q € D; 



(1.5) 



(1.6) 



that is, M — Mi. Here, each Zi is a positive constant chosen so that Mi — 1 (we can 

do so because of Hypothesis IaI. below). Thereby, /q M = 1. We note that since Ui is assumed 
to tend to -l-oo as approaches dDi, the corresponding partial Maxwellian Mi tends to as 
approaches dDi, i e [N]; consequently, M tends to as q approaches dD. The fact that the 
Maxwellian factorizes — which comes from the fact that the energy stored in the chain is the sum 
of the potential energies stored in each spring — will be crucial throughout the rest of this paper. 
For a start, this fact allows us to write 



(1.7) 



Multiplying (ll.4al) by f/M, using (II. 7p and (formally) integrating by parts, the corresponding 
weak form of (|1.4p is: Find -0 = "0(3;, q, t) such that 

V i — 1 j — 1 

4Wi(iV-fl) ^ I ^ ' 



for all ip — (p{x, q) in a suitable function space. 
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For the sake of convenience we define the foUowing bilinear forms: 

r(«;.,.) -X^/iv.(-)^, ^.r) - i^X^^V... (1.9) 



Tiu; a, t):^- f ^(V,M)q,a • V,^ (^) , (1.10) 
^(-)-I„Ei:i|MV,©.V„(^). (l.U) 



Then, (II. 8p can be written concisely as 
'9V 



, ^/M^ + T{u; ip, ^) + /C(^, ^) + r(M; <y5) + /C(7A, (^) = (1.12) 

for all iy9 = V?(ic, g) in a suitable function space. We note that T and JC involve partial derivatives 
of their arguments with respect to the spatial variable x only. Analogously, T and /C involve 
partial derivatives of their arguments with respect to the configuration space variable q only. This 
motivates the use of the alternating direction scheme based on operator splitting whose informal 
description is given in the next subsection. 

1.3. Alternating direction scheme. Let At be such that M := T/Ai e N and define :— nAt 
for n € {0, . . . , M}. We will consider the following alternating- direction semidiscretization of 
(|1.8p : We initialize the scheme by defining ■— ''Po', for n g {0, . . . , A/ — 1} and then define the 
'intermediate' function and the approximation ^/;"+^ to ■!/'(t"+^, •, •), respectively, by 

( "^"X^"^" ' M ) + ^) + ^(^"^'/^ ^) 

= -r(w(-, r); V^", ^) - IC{r,f) (1.13a) 

and 

/ n+l _ „;.n+l/2 \ 

,S)+/c(v/'+\^)--r(«(-,n;V/\^) 



At/2 ' M 



f{u{., t"+i); ^) - ^(^"+1/2, ^), (1.13b) 



for all (p = (p{x, q) in a suitable function space. In (jl.l3ap the spatial bilinear forms T and /C are 
treated implicitly while the configuration space bilinear forms T and /C are treated explicitly. In 
(I1.13b|) the spatial bilinear forms T and /C and the configuration space bilinear form T associated 
with the drag term are treated explicitly, while the bilinear form K, is treated implicitly. 

Let ( {q^''\ w'^'^ )j and ( {x^'^\ w^^'' n be ^j^- and 1-weighted quadrature rules on D 

V / kelQo] ^ '' ke[Qs i] 

and Q, respectively. We then approximate (|1.13ap by performing numerical integration over the 
configuration space, which results in 



fe=i 



+ J2wl,'^ I div. («(, r+i)^"+i/2(., )) ^(., gW ) 
fc=i -^^^ 

Qd i- N 

5:4^) / EM('7^'^-^)(V.n(-,n)gf ^"(•,g«)-V,^ Q 

7 1 ^ r,' 1 
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k=l 



N N 



n 4Wi 

" i=i i=i 



V/' 



• V 



(..qC^)) 



(m) 



for all — ip{x, q) in a suitable function space. Here, the symbol denotes equality, up to 
quadrature errors. By selecting Qd linearly independent functions C(m)j G [Qd], of q G D 
such that C(m)(9*''^'') = ^fcm, k,m £ [Qd], and taking successively c/? = t/?(m); where (a;, q) := 
x(^)C(m)(9)j in thfi equality above, we obtain a total of Qd independent variational problems, 
each posed over the c?-dimensional domain f2, of the form: 



A</2 



4Wi(iV- 



-, Qd r /• ^ 

Wd fc=l L-^" ^=l 



C(ji 



• V 



1i 



C(m) 

M 



M 



(■,'?('=>) 



=:OT(™)(^";x) Vme[QD], (1.14) 



for all X — xi^) in a suitable function space, where each 9Jl(m)('0"; •)' ^ [Qd], is a linear 
functional. Thus, (jl.l4p amounts to solving Qd mutually independent linear convection-diffusion 
problems over fl. 

In turn, we can approximate (jl.l3bp by performing numerical quadrature over $7, resulting in 



E 

fc=i 



w, 



(k) 



At/2 



M 



/C=l "^^ 1=1 7 = 1 



E^ 

k=l 



(fc) {Iq/LqY 



4Wi(iV+l) 



n+l/2 



1 

(x(*),.) 77' 



for all (y9 = q) in a suitable function space. By selecting Qji linearly independent functions 
X(jn)^ m e [Qn], of a; S 17 such that X(jn) (^j'"'^'') = ^km, k,m & [Qn], and taking successively ip — 
(p(jn), where Lp(^„^){x,q) := X(m)C('7), in the equality above, we obtain a total of Qn independent 
variational problems over the A^d-dimensional domain D of the form: 



M 



Vq, ( 4 



M 



- /■ div. («(.,r+i)^"+V2^ 



w, 
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(fc) (Iq/Lq? f n+1/2 . V ^ 



I (_ 



=:9T(™)(7A"+i/^C) ymeiQn], (1.15) 

for all C = ((q) in a suitable function space, where each Vi(^rn)i'^^~^^^^', G [Qn], is a linear 

functional. Thus, (|1.15p amounts to solving [Qn] mutually independent linear elliptic variational 
problems, each posed on the high-dimensional configurational domain D = Di x • • • x Dn C M.^'^. 
It is the approximate solution of (|1.15l) by greedy algorithms that this paper is concerned with. 

2. The configuration space operator 

2.1. Variational formulation and function spaces. The form of the problem (|1.15|) motivates 
us to consider the linear elliptic variational problem 

ai^,^)^fiip), (2.1) 

posed on the high-dimensional configurational domain D = I?iX---xL)jvC M^'^, where 

/„|:i:3l;Mv„(^)-v„©.e/^l^, p., 

i — 1 j — 1 

the parameter c is positive and / is a linear functional. The natural function space associated 
with problem (|2.ip is 

H(D; M) e Ll/^iD) n M LL(D): Vq^(^/M) e KiD)f Vz G [TV]}, 

equipped with the norm 

N . 1/2 
|2 ^ 



II¥'IIh(D;M) := ( ll^llL?,JD)+Ell^'J.(^/M)||[L^(Dr 
^ i=l 

The spaces L^y^(D) and H(D; M) are isometrically isomorphic to, respectively, L^(D) and 
Hjyi(D) via the relations 

l2/^(D) = ML2^(D), |1.|1j^,^^(^^ = ||M-i.||l,^^^, (2.3a) 

H(D;M) = MHi,(D), IHIh(D;M) = ||M-^-||h1(d)- (2-3b) 

Later, we will make use of the spaces ll{Di; Mi), i E [N], each of which is the i-th partial 
Maxwellian analogue of H(D; M). That is, 

H(A; Af,) := e L2/^,,(A) n M,LL(A): Viip/M,) G [Ll^iD,)]^], 

I 1 2 

equipped with the norm ||</?||H(D,;Af,) — \ II'/'IIl^ (d,) + ^^l^'^^)\\\\.l,\D,)Y 



Remark 2. 

(1) For i e [N], R{Df,M,) is exactly H(D; M) if iV = 1 and M = M,. None of the results 
involving H(D; M) appearing below depend on restrictions on N and thereby remain valid 
for ii{Di;Mi). Just like (j2.3p . (p i— Mi(p is an isometric isomorphism between L\,j,{Di) 
and Lyj^j,{Di) and between ll\j,{Di) and H(_Di; M^). 

(2) The definitions above can be extended to open subsets of D and of the Di, i e [N], in the 
usual way. 

Before listing our structural hypotheses and proving the properties we need of H(D; M) we fully 
state the weak formulation of our model problem: 

Given / e H(D; M)', find V e H(D; M) such that 

ai^P,ip)=f{cp) V^eH(D;M). (2.4) 

We adopt the following structural hypotheses. 
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Hypothesis A. For each i e [N], the sprmg potential Ui belongs to C^([0, y)), where bi > 0, 
and satisfies lims_j.f,./2_ U{s) — +00. 

Immediate consequences of Hypothesis IaI are that M € C(D)nC^(D) and that, for any if (S D, 
there exist positive constants ck and Ck such that ck < M(qr) < Ck, for all q E K. 

Hypothesis B. For each i £ [N], H]^,^. (i?i) is compactly embedded in Lj^,{Di). 



Remark 3. It is easy to check that springs obeying any of the example force models (|1.2p and (11.31) 
comply with Hypothesis 0- 

In Step 1 of section A.l of 'BS08| it is proved that springs obeying the FENE model (|1.2I) satisfy 
Hypothesis IbI. under the condition hi > 2. The compliance with Hypothesis [bI of springs obeying 
the CPAIL model p.3p is shown in Lemma A.l in Appendix under the condition hi > 3. 

Lemma 2.1. L^(D), HjJf(D) for m £N, Lf^j^{D) and H(D; M) are separable Hilbert spaces. 

Proof. The operation Lp e L^(D) i-> Lp/y/M defines an isometric isomorphism between L^(D) 
and L^(D). Therefore the first space inherits its separability from the latter. On noting that 
G Lj'Qj.(D), Theorem 1.11 of |K084) guarantees the completeness of H}Jf(D) (this source 
actually states the result for the case m = 1 only; however, the proof carries over to higher m in 
this single- weight case) and thus, HjJj'(D) is separable by an argument along the lines of |AF03[ 
^3.5]. The spaces L^^|^(D) and H(D; M) inherit these properties via the isometric isomorphism 
(|2.3p . Finally, as their respective norms obey the parallelogram law, these spaces are Hilbert 
spaces. □ 

Lemma 2.2. H^(D) is compactly embedded in L^(D), and H(D; M) is compactly embedded in 
L?/m(D). 

Proof. Throughout this proof we will assume, for ease of exposition, that N — 2] the argument 
carries over to higher N without difficulties. Let u £ H^(D). As, by pT^ . M = Mi (g) M2, it 
follows from Fubini's theorem that, for almost all £ Di, 

u{qi, •) £ Lm,{D2) n Ll^{D2) and <9aM(gi, •) £ L^,^^(L»2), 

where a is any multi-index in [No]'' with < |a| < 1. Fubini's theorem, again, ensures that, given 
ip2 £ Cg°(i?2) and a2 £ [Nq]'^, < aa < 1, 









J D 



d(0,a2)u{qi,-)V2 



(^idqi. 



for all ifi £ Cg°(Di). Therefore, da^iuiq^ •)] = c^(o,q2)^(9i' ') weak sense on D2 for almost 

all <7i £ Di. As 9(o,Q2)'"(<7i, •) lies in 1j\j^{D2) for almost all q^ £ Di, we have that 

u{qi, •) £ ilM2{D2) for almost all q^ £ Di. (2.5) 

In the same way it can be proved that u{-, (72) £ H]^,J^ (Di) for almost all 92 G ^2- 

Let us define, for i £ {1,2}, the sequence (Di^(„))^^^ of bounded and proper subsets of Di by 

A,(n) :=i?(0,^). Then, 

00 

C A,(n+i), n£N, IJ = Di and Hj,^. (i?j^(„)) d L^^^ (i?j^(„)). 

Tl=l 

This last relation is a consequence of the corresponding relation for the unweighted case, H^(£'i_(„)) 
(S L^(Dj („)) — in turn a consequence of the boundedness and Lipschitz continuity of i?i,(n) — on 
account of the existence of positive lower and upper bounds for Mi on Dij^^)! whereupon there is 
algebraic and topological equivalence between ^\[. (Dij^n)) and H^(i3j („)) and between L^^. (i?i^(n)) 
and L2(D,_(„)). 

Letting, for n £ N, D(„) :— X^^^ £ the above properties get inherited: 

CO 

D(„) C D(„+i), n e N, y D(„) = D and H^D(„)) d L^(D(„)). 

n=l 
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The third statement follows from the fact that the D(„), being Cartesian products of bounded 
Lipschitz domains, are also bounded Lipschitz domainfQ. Let us define Di \ Dij^n) and 

D(") := D \ D(„). Thanks to |OK90I Theo rem 17.6], the above compact embcddings on members 
of a nested covering imply the following characterizations (the first, for i e {1, 2}): 



(A) Ll,^ (A) ^ Hm sup / uHd, I \\u\\L (^^^ 
Hm(D) (sL2,(D) ^ hm sup / w^M/ IIt^II^wd) =0- 



= 0, 



(2.6) 
(2.7) 



From Hypothesis 0, the left-hand side of (|2.6p holds; hence, its right-hand side also holds. Using 
(|2.5p and (|2.6I) with i = 2, we deduce that for each e > there exists some h = n(e) S N such 
that n > h implies 



/ u^M= [ [ 

JDixD<"> Jdi Jl 



Mi(q-i)dqi 



l|w(gi,-)llHi (D2)^'^T-iii)^ii 

JDi ' 



1 UiQir) 


+ |Vq^u(qi,-)|') M2 


IJD2 ^ 





Afi(qi)dqi < e|lu|lH^(D) . 

An analogous result can be proved for the M-weighted integral of v? on D^l'^ X D2. Then, since 
D(") = {Di X Z)^"^) U (L»J"^ X Ds), the right-hand side of ^TT^ holds; hence, so does its left-hand 
side. 

Finally, the embedding H(D; M) d L^y|^(D) follows directly from the embedding H||^(D) <e 
L^(D) on account of the isometric isomorphism (|2.3p . □ 



Lemma 2.3. The following inclusion holds: Co(D) C H(D; M). 

Proof. Let ip G Co(D) and K :— supp((^) d D. Then, trivially, </? G L^^,y(D), since 



M 



^2l<|if|sup^<oo, 
A- M - ' M(9) 









M 



which, in turn, stems from the fact that M is positively bounded from below on each compact 
subset of D. Similarly, for all iiT' d D, 

<\K'nK\ sup M£^<oo 
qeK'nK M(q) 

on account of which <p g The latter implies that ip/M defines a regular distribution 

in the usual way. Then, for each i G [iV], Vq. ((^c/M) exists as a distribution and coincides with 
the classical «-th component gradient of (ys/M, which belongs to [C(D)]'' because of Hypothesis IaI. 
Then, 



M < \K\ sup 

gS-ff 



1(g) 



M(q) < 00 



and that proves the lemma. 



□ 



This follows by combining Theorem 3.1 in the Ph.D. Thesis of Reinhard Hochmuth: Randwertproblem einer 
nicht hypoelliptischen linearen partiellen Differentialgleichung, Dissertation, Preie Universitat Berlin, 1989, which 
implies that the Cartesian product of a finite number of bounded domains, each satisfying the uniform cone property, 
is a bounded domain satisfying the uniform cone property, and Theorem 1.2.2.2 in the book of Grisvard IGri85| . 
which states that a bounded open set in R" has the uniform cone property if, and only if, its boundary is Lipschitz. 
In the special case of the domain D(„) an alternative proof is to note that, as a Cartesian product of bounded open 
convex sets, D(„) is a bounded open convex set in R" (cf. IHULOl| . p. 23), and then apply Corollary 1.2.2.3 in 
Grisvard IGri85l , which states that a bounded open convex set in has Lipschitz boundary. 
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2.2. Properties of tensor products. 

Lemma 2.4. Suppose that T G 2?'(D) is a distribution such that 



N 



T (^0 (^('' ) = 



ie[N\ 



Then, T ^0 in V'{D). 

Further, for any ensemble of sequences of distributions [Rn^j , i € [N], with Rn^ G V^Di 

V / n>l 

and such that lini„_j.oo Rn^ — -R*"*' in T)' [Di) for i G [N], we have that 



lim 



Ue[Af] 



(g) 

ie[N] 



m V'{D). 



Proof. These are standard results from the theory of distributions, so we omit the proofs and refer 
the reader to Section 1.3.2 of the book of Vladimirov |Vla02j . for example. □ 



Lemma 2.5. The following statements hold: 

(1) For any ensemble r^*' G }i{Di;Mi), i G [N], (g) 



(2) Suppose that r^^ : A - 
ments are equivalent: 

(a) r(^) G H(A; M^) \ {0} for all i G [N] 



ie[N] 



gH(D;M). 



, i G [N], are measurable functions. Then, the next two state- 



(b) (g),, 



e[N] 



gH(D;M)\{0}. 



Proof. (1) It is immediate from the factorization of M that {S)jg[jv] 
to Lemma [23, the identity 



M 



N 



m7 



belongs to (D). Thanks 



(2.8) 



holds in the distributional sense. Then, as r^*' /M^ G L|,^. (D^) for i G [A^] \ {j}, and V{r'^^^ /Mj) G 



[hl^XDj)]'', the factorization of the Maxwellian M implies that, for j G [A^], Vq^.(0,^^ r'-''> /M) G 
[L^ (D)]'^. That completes the proof of Part (1). 

(2) We shall prove the second part by showing that (b) is both necessary and sufficient for (a), 
(a) => (b): This is immediate from the first part and the fact that the tensor product of the 



rii) 



, i G [N], cannot be null if none of its factors is. 
(b) = ■ 



(a): Suppose that 0iLi t"^'-* G H(D; M) \ {0}; then, because of the tensor-product 
structure of M, the positivity of Mi on compact subsets of Di for i G [N] and Fubini's theorem, 
r(*) G AhLl^iD,) nLl^{D^), i e [N]. Hence, each r^^V^i defines a regular distribution in 



T>'{Di). Again, Lemma l2.4l makes (j2.8p valid and thus. 



N 


2 N 






i=l 


H(D;M) 



N 

E 



N 

n 



j = l li=l 



(2.9) 



Now, none of the r*^*^ can be null (otherwise their tensor product would be null). On combining this 
with their l/M^-weighted square integrability, the identity (j2.9p yields || V(r(*YAfi)|| [l2^ (-Di)]'' < °° 
for all i G [iV]. Hence r^^ G H(A; A^j) \ {0} for i £ [N]. ' □ 



3. Separated representation 

3.1. Two algorithms. The existence of a unique weak solution to (|2.ip is an immediate conse- 
quence of the Lax-Milgram theorem via the facts that H(D; M) is a Hilbert space (cf. Lemma |2. 11 ) 
and a is a bounded and coercive bilinear form on H(D; M). By virtue of the Riesz represen- 
tation theorem, there exists a bounded linear operator A: H(D; M) — ^ H(D; M)', defined by 
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{Aip){ip) = a{ip,(p) for all ip E H(D; M). Thanks to the symmetry of a, the weak formulation 
(|2.ip can be restated as the following, equivalent, energy minimization problem: 

:= argmin J/((/5) where J/ ((^) := ^ a(<^, (/j) - /(</?). (3.1) 

V5eH(D;M) ^ 

We observe that, with if; g H(D; M) as in ([XT]) . 

Jf(^)^laicp-^,cp-i:)-^ai^,iJ;) V^eH(D;M). (3.2) 

Following the work of Le Bris, Lelievre and Maday |LBLM09] concerning the numerical solution 
of high-dimensional Poisson equations, we consider two numerical methods. 

Algorithm 1. (Pure Greedy Algorithm) 

0. Define: fo := / € H(D; M)'. 

1. For n > 1 do: 

1.1 Find rli^ G H(A;Mj), i £ [N], such that 

/ ^ 

r.(i) 



(r^l^...,rrOe argmin J/„_, • (3.3) 



1=1 



(s(i),...,s("))eXf^i H(Di;A/,) 

1.2 i^e/ine.- /„ := /„_i - ^ ({gj^^ ri^^) e H(D; M)'. 

1.3 // ||/n|lH(D M)' — TOL, then proceed to iteration n + 1; else, stop. 

Algorithm 2. (Orthogonal Greedy Algorithm) 

0. Define: /o := / G H(D; M)'. 

1. For n > 1 do: 

1.1 Find r'h^ G }i{D,;M,), i £ [N], such that 

(rW,...,rW)G argmin J/„_/(8) r«) . (3.4) 

(s(l),...,s<"))eX" 1 H(Di;J\/i) \«=i / 

1.2 Minimize J r on the span of r^i''] : i.e., find a'-"-' G R" such that 

V /fee[n] 



a 



= arg min J/ ( ^ r« ) . (3.5) 



/3GK" i=l 



1.3 Define: U -.^ f - A (ELi 4"^ 0^1 rf) G H(D; M)'. 

1.4 // ||/n|lH(D-M)' — TOL, then proceed to iteration n + 1; else, stop. 

For future reference, we define ipn G H(D; M) as the unique solution of the problem 

a(^„,^) =/„((p) V^GH(D;M). 
Clearly, for all n up to the (existing or not) termination of the corresponding algorithm. 



V'n-i — ^iLi foi" the Pure Greedy Algorithm, 

ip — J2k=i Q^fe"^ 0iXi ^fc ^ fo'" the Orthogonal Greedy Algorithm, 



where ip = ipo is the unique solution of (13. ip . Proving the convergence of the algorithms amounts 
to showing that the sequences {tpn)„yQ defined by p.6p converge to in H(D; M). 
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3.2. Correctness of the algorithms. The proof of the correctness of Algorithm [l| (respectively 
Algorithm Q) amounts to showing that, given /„_i € H(D; M)' (respectively (/„_i, a'"^^)) G 
H(D; M)' X R"-i), the loop 1 returns a well-defined member of H(D; M)' (resp. H(D; M)' x M"). 

We start by observing that, thanks to the first part of Lemma I2.5L the set of A^-way tensor 
products of ensembles of functions H(Di; Mi), i S [A^], is a subset of H(D; M), thereby rendering the 

minimization problems (13. 3|) and p. 41) sound. However, the existence of sol utio ns {rn^^ , ■ • ■ , ) 
to these problems is quite another matter: it will be proved using Lemma l3.ll and Theorem |3.2 
below. 

Lemma 3.1. Suppose that f £ H(D; M)' \ {0} and consider the functional Jf, as in (j3.ip . Then, 
there exists (r(^\ . . . , r^^^) in xf^^ll{Di; Mi) such that 

Jf (^(8)r«^ <0. 

Proof. Consider any functional / £ H(D;M)' \ {0} and assume that the thesis is false; i.e., 
J/(0^irW) >Ofor ah ensembles (r^^), . . . , r(^)) e X H(A; M,); then, 



1 

2" 



N N \ / N \ N 

0r«,0r« >/ {8)r« V (r^ , . . . , r^) G Xil{D^■,M,). 

i=l i=l J \i=\ / i=l 



Given a particular ensemble (r*^^^ , . . . , r'^^^) e X H(£'i; Mi), we can replace r'-^-' with er*^^-' and, 
by virtue of the bilinearity of a and the linearity of /, we obtain 

1 / ^ ^ \ / ^ \ 

-e^a 0r-«,0rW >e/ 0r« . (3.7) 

^ \i=l i=\ J \i=\ / 



By combining the inequalities resulting from dividing both sides of (j3.7p by positive £ and taking 
the one-sided limit e — ?> 0-(- and from dividing (j3.7l) by a negative e and taking the one-sided limit 
e — >■ 0_ we get that 

/(^(|r«) =0. 

As this is valid for any ensemble (r^^-*, . . . , r*^^^) e X^j^ H(Z3i; Mi), Lemma 1 2. 3l impli es th at it is 
valid, in particular, for any ensemble (r(^\ . . . , r*^^^) e X^j^ Cg°(Z?i), whence Lemma 1 2. 4 implies 
that / = 0. As this contradicts the hypotheses of the lemma, its thesis holds. □ 



We are now in a position to prove the existence of solutions to problems (|3.3I) and (j3.4l) . 

Theorem 3.2. Given fn-i £ H(D; M)', each of the problems (|3.3I) and p.4p /las a solution. 

Proof. Since problems p.3p and p.4[) are completely analogous, it suffices to consider one of 
them — say, p. 31) . Then, as (0, . . . , 0) is a solution of p.3p and (13.41) when /„_i = 0, we assume 
from now on that fn-i 7^ 0. 

By p. 21) and the coerciveness of a, Jf„_i{(p) > ~^a{tp,ip) for a ll ip £ H(D; M), where ip is the 
unique solution of (|2.4I) in H(D; M) when / = /n-i. As, by Lemma |2.5| . the A^-way tensor product 
of functions in H(Di; Mi), i £ [N], is a subset of H(D; M), J/„_i is bounded from below over that 
manifold. That is, 

m := inf J/„_, ((g) s(*M > -oo. (3.8) 

(s(i),...,s("))eX" 1 H(Di;Mi) Vi=l / 



It follows from Lemma l3.ll that m < 0. Our aim is to show that the infimum m is attained at an 
element of the form rW with (r^^), . . . ,rW) £ xfL^{}i{Df,M,) \ {0}). 

From p.8p , there exists a sequence (^(g^g [^v] ^i*^ ) N-way tensor products of functions in 

H(A; Mi), i£ [N], such that 

lim Jj^_, (^r^A =m. 
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On noting that, from the definition of a in (j2.2p . for all (p E H(D; M), 



1 . 

> — mm 

- 4 



An 



4Wi 



H(D;M) 



a(-0,V'), 



it follows, by setting — <^ 



ie[N\ """k ' 



that the sequence ( 



>>ie[N\ ''"k 



(0 



k>l 



is bounded in H(D; M); 



in other words, there exists C > such that (cf. (|2.9 



i=l 



H(D:M) 



N 

n| 

i=l 



N 



N 



E ni 



[LL.(i).)l'' 



< C 



(3.9) 

for all fc > 1 . Since the value of rj, is unaltered by multiplying the first iV — 1 factors 

by positive constants ci^^, . . . , CN-i,k, respectively, and dividing the final factor by the product 
ci,fc • • • CN-i.k, we can assume without loss of generality that 



ADi) 



1, ic[N^l]. 



(3.10) 



Thus, it follows from that 



IN) 



IN) 



N-1 



[LL,(I).)l'' 



V(rfVMA^) 



[LL„(Civ)]'' 



< c. 



(3.11) 

is bounded in H(D; M), and H(D; M) is a Hilbert space, 
and therefore reflexive, the sequence has a weakly convergent subsequence in H(D; M), denoted by 



Since the sequence ^0, 



ie[N] ' k 



k>l 



(0. 



''ie[N] ' (i>{k)J ' denote its weak limit by r € H(D; M). Since J/„_i is convex on H(D; M) 

and continuous (and thereby also semicontinuous) in the strong topology of H(D; M), it is weakly 
lower-semicontinuous on H(D; M). Hence 



J/„-i(r) < liminf J/„_, 



fc— f C30 



N 



N 



= lim J/„_, ( (g) r^'M = m < 



Thus we deduce that r ^ (as r = would imply that Jf„_i{r) = 0); hence, r E H(D; M) \ {0}. 
According to (|3.10l) and (13.111) each subsequence i^^^^k) ^ _ j is bounded in the respective space 

lj\^^_jiDi), for i e [N]. Then, (^'"^(fe)) has a weakly convergent subsequence in L^^j^^. (D^), say 

(^^'''(fc)) ' '^'-'^ * ^ I^]' -^^^ denote by r^'^ e L^^j^^ (1?^) the corresponding weak limits. As by 
LemmalM G^{D^) C H(A; M,) C L2/^,^^ (A), for aU ^ G Cg°(A) the mapping ^ e L2/j^,^^(A) ^ 
(v'l OlJ^jj- (_Di) defines a bounded linear functional on L^^^,^ (Di). Thus, {j'^^{k)l^-^ij _ converges 



to r(*) /Mi in 2?'(£)i) for i e [TV]. Hence, by Lemma |2.4 



N r 



0'(fe) 



in 2?'(D). 



(3.12) 



Similarly, the inclusion C3°(D) C H(D; M) (cf. Lemmal2J) and the fact that, for all ip e C[f (D) 
the mapping ^ £ H(D; M) t-^ {f,OLj^^{D) defines a bounded linear functional on H(D; M) imply 



lim (8)4^ 



lim 



0ti 



N ^(*) 



0'(fc) 



in ^'(D) 



(3.13) 
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on account of r being the H(D; M)-weak limit of the sequence (^{^j^j^v] 4>(k)) ■ 'D'{D) is a 
Hausdorff topological space, the limits in (j3.12l) and p.lSp have to coincide. That is, 

N 

M- V = (g) r-W inX''(D). 

Hence, r = almo st ev erywhere. As r e H(D; M)\{0} and has a tensor-product structure, 

the second part of Lemma implies that r^'^ e H^Df, Mi) \ {0} for i e [N]. Now, 

Recalling the definition of m from p.Sp . we have thus shown that the infimum in p. 81) is attained 
at {g),^irW. Thus, (r^^), . . . , r^) e X^i(H(A; M,) \ {0}) is a solution to problem 

□ 

Having proved that the minimization problems (|3.3p of Algorithm [l| and (13.41) of Algorithm 
have solutions, establishing the correctness of what is left of the algorithms is straightforward. 
The Galerkin problem 1 . 2 of Algorithm [i] is well-defined and has a unique solution for each n > 1, 
because it is equivalent to the minimization of a coercive quadratic form over a finite-dimensional 
linear space. Then, at last, the definition of the n-th residual in step 1.2 of Algorithm [l| and in 
step 1 . 3 of Algorithm [2I are correct on noting that A maps H(D; M) into H(D; M)'. 

In the next section we establish the convergence of the two algorithms. 



4. Convergence of the Algorithms 
4.1. Euler Lagrange equations. 

Lemma 4.1. Local minimizers {rn \ . . . ,ri^'') of the minimization problems (j3.3p or p.4p satisfy 
the following Euler-Lagrange equation system: For all {s'^^\ . . . , s^^^) € y<ie[N] H(I?i; Mi), 



From this, it follows that, for the Pure Greedy Algorithm (Algorithm]^): 

a(^n,&rl:A =0. (4.2) 



Proof. Let {rn \ ■ ■ ■ , r'i^^) be a solution to the minimization problem (j3.3p or (j3.4p . Then, given 
any ensemble (s*-^-*, . . . , s^^-*), (|4.ip is but a way of writing that the derivative of 



N , 



with respect to e is zero when evaluated at e = 0. As, by hypothesis, (r^n' , ■ ■ ■ ,r^n^^) is a local 
minimizer of J/„_i and s 1— >■ 3n(e) ■— Jf„-i (^(S'iLi {j'n^ + es*^'-'^^ is regular enough, the fact that 
j;^(0) = implies that (gD) holds. 

Setting (s*-^-*, . . . , s^^-') — {rn \ . . . , ri^-*) in (14. ip . we obtain from the definition of the i/'n that 

N N \ / N 



a (^grW, rWj ^ a (^^„-i,grWJ . (4.3) 
Combining this with (j3.6p we obtain (|4.2p . □ 
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Remark 4. 

(1) The above lemma only states that local minima of the minimization problem p.3p and 
(j3.4p satisfy the Euler-Lagrange equation (|4.ip . The converse may be false, of course: 
although the functional that is minimized is quadratic, the set over which it is minimized 
is nonlinear, so there is no reason why a stationary point should be a local minimum. 

(2) In what follows we make liberal use of the norm ||'||^ := a(',')^/^ on H(D;M), which, 
thanks to its equivalence with |HIh(d m)' makes no difference when making topological 
statements (such as convergence). 

Lemma 4.2. Let (rn^^ , ■ ■ ■ , J'i^'') be a global minimizer for the minimization problem 
Algorithm]^. Then, 



N 

(2) 



(V'n-l,0tl^" 



Proof. The first equality in 
be written as 



comes directly from 



sup 



a(V'„-i,s) 



Now, analogously to 



of the 

(4.4) 
Jfn~i can 



JU-lif) = ^O-i'f - ^n-l,<y5 - ^n-l) " ^a(V'„-i , "^n- 1 ) V G H(D; M). 

Combining this representation of with the fact that r„ := r„ minimizes J/„_i among 

the members of {S)ie[Ar] H(Di; M^) and the first equality of (|4.4p . according to which a (7/)„_i, r„) = 
||r„||^, we have, for ah s e 0,grjvi H(A; M^) \ {0}, that 



elN] ■ 
a ('(/'«- i,r-„) 



< 



a(V'„-i,s) 

Vn-l —75 S 



Therefore, 



a{rn,rn) a(s, s) 

Taking the supremum over s G {S)ig[jv] H(-Di; Mi) \ {0} and noting that r„ is an admissible s we 
get the second equality in (14.41) . □ 



4.2. Convergence. 

Theorem 4.3. T/ie Pure Greedy Algorithm (Algorithm]^) converges to the solution tp to (j2.4p . 



Proof Let (( 



„(i) 



be a sequence in X^2H(£>i;Mi) returned by the Pure Greedy 



Algorithm and let us adopt the shorthand notation r„ ;= 
in Lemma 4.1 we obtain 



Then, from ([H^l) and (02]) 



ll^n-llla = W^n + rnWl = \\^n\\l + IKWl ■ 



Hence the sequence (||'*/'ti||£i)„>o is nonnegative and monotonic nonincreasing, and therefore con- 
verges in K; by summing the above expression over n we then deduce that 



'^a{rn,rn) < 00. 



(4.5) 



Let us define the function 0: N ^ N recursively by 4>{1) := 1 and 

(j){k) := argmin {||r„||^ < Hj'^Cfc-i) D , k>2. 

n>ct>(k-l) 

From (|4.5p the function (f) is well-defined and strictly monotonic increasing. Hence, it is suitable 



for defining subsequences. As each (' 



,(1) 



(N) 



) is a global solution to the problem (j3.3l) with 
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the instance fn-i, via p.6p and Lemma l4 .21 we have, for n>m> 1, 



||V'0(n)-l - '00(m)-l||„ = ||V'0(n)-l||(j + 1 1 V-'^i(m)- 1 1 1 - 2a I V'0(n)- 1 , ''/'^(n)- 1 + X/ 

\ k—(p{m) 
0(n)-l 

= ||^0(m)-ilL- llv-^M-ilL-z J2 «(V'</.(«)-i,f-fc) 

0(n)-l 

< IIV'^M-ill' - llv-^cnj-ill' + 2 ^ IkfclLlk^wlL 

0(n)-l 

< IIV-^lnO-ilL - IIV-^H-ilL + 2 ^ llrfcll^. 

fc— (/)(m) 

From the convergence of (||^0(n)-i in ^ ^nd (|4.5p . we deduce that the sequence ('00(n)-i)„>]^ 

is a Cauchy sequence in H(D; M) and thus converges to some ipao & H(D; M). Another consequence 
of the global optimahty of each (ri^\ . . . , ri^^) is: For all (s^^^ . . . , s^^^) G X ,;g[7v] H(Di; M^) and 
n>l, 



1 / TV N 



i=l i=l 



a (^^4>{n)-l, > Jf4>(^)-1 (^0(n)) 



-a (r0(„),r0(„)) - (^^(n)) 
-^a(»-0(n),r^(«)) ■ 



Taking the limit as n tends to infinity at both ends, and noting that by (|4.5p the right-hand side 
of the last inequality converges to 0, we obtain 



1 / AT AT 



a(^^oc,^y''>^ >0. 



Thus, Lemma 3A implies that ipao = 0. Hence the sequence (||''/'0(n)-i converges to zero 
as n — > oo. As the sequence (||V'n||£i)„>o is monotonic nonincreasing and {<j){n) — l)„>i is a 
monotonic increasing infinite sequence in N, if follows that the full sequence (||V'n||)„>i converges 
to the common limit in R: = HV'ooHai giving lim,n.oo V'n = in H(D; M). □ 



The following corollary is a direct consequence of Theorem 14.31 and will prove useful later on. 

Corollary 4.4. Let Fi he a dense subset ofll{Di;Mi) fori £ [N]. Then, the span of^^^^j^-^Fi 
is dense in H(D; M). 



Proof. Let r G H(D; M). Applying Theorem |4.3| to the case in which the right-hand side functional 
/ e H(D; M)' of problem (12.41) is H> a(r, ip) (i.e., the H(D; M) approximation problem) it follows 
that T can be approximated arbitrarily closely by finite sums of the form X^melA/l^j "^^^^ 



where M e N and rl^ € H(A; Mi) for m e [M] and i e [N]. Thus, if we can show that 0igrjvi Fi 



e [M] yyie [n] ' ™ ' 
e[N] 

is dense in the manifold 0jg[jY] K{Di;Mi), our desired result will stand. 

Let, then, r(') G H(A; M,), for i e [N]. From the density of F, in H(A; Mj) for each i € [N], 
there exists a sequence (rn^j in Fi, which converges to r'^*^ in H(Di;Mi). Now, 

V /n>l 



N N 

{g)rW_0r«=^0t«, where t« := <j _ if , = fc, 



r„ 



if j > fc, 
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Then, (cf. (EH)), 

II'^"IIh(D;M) — ^ 



N 



k=l 



N 



N N 
j=l i=l 



As each product term on the right-hand side above consists of — 1 bounded factors and one 
vanishing factor as n — > oo, the fuh expression tends to zero as n tends to infinity and, therefore, 
so does the left-hand side. The desired resuh follows. □ 

Remark 5. Suppose that, for each i G [N], 

Cg°(A) is dense in H(A;Mi). (4.6) 
Then, as span Cg°(A)) C Cjf (D) C H(D; M), we have, thanks to CoroUaryO, that 

CS°(D) isdenseinH(D;M). (4.7) 

Springs obeying the FENE model (jl.2|) comply with (14. 6p under the condition &i > 2 as is proved 
in Remark 3.7 of I'MasOSI. Springs obeying the CPAIL model (|1.3p . in turn, comply with (|4.6p as 
it is shown in Lemma I A. ll in Appendix I Al. under the condition bi > 3. So, in these two cases, (14.71) 
holds. 

Interesting as (j4.7l) is, we make no use of it in this work and that is why we shall not adopt 
(|4.6|) as a hypothesis on a par with hypotheses \X\ and iBl above or hypotheses ICl iDl and [El below. 
However, we do use (|4.6|) as an ingredient of the proof of the compliance of FENE and CPAIL 
spring potentials with Hypothesis [c| of Section [HI (cf. Corollary |C.2| in Appendix 0) • 

Theorem 4.5. The Orthogonal Greedy Algorithm (Algorithmic converges to the solution ip to 
problem (|2.4|) . 



Proof. We first note that thanks to (j3.6p . the optimality of in (1X3 and the optimality of 



\i n T ■ ■ ■ , I n 



) in (|3.4p (via Lemma l4. II). 



fc=i 



' k 

= 1 



< 



N 



l'0n- 



N 

i=l 



Thus, just like in the proof of Theorem l4.3l . we have that the real sequence (||'0n||Q)„>o decreasing 



and thus convergent and that J2n>i ^ \ ® 



[AT] T-n ;0ig[Ar]''n j < OO. As (V'n)„>o ^ bounded 
sequence in the Hilbert space H(D; M), a weakly convergent subsequence (V'0(n))„>]^ can be ex- 

_i_j^) with 



tracted and we denote the weak limit by tpoo- From the optimality of (' 



,(1) 



(n)-l-l' 



respect to problem p. 41) we have by Lemma l4.ll that, for all {s^^\ . . . , s^^^) G X j^j^v] H(Di; Mi), 



1 



N 



N 



(0 _ 



vi=l i=i 

Taking the limit ti -> oo at both sides yields 



a (^^^(„), > --a (^grW„)+i,gr«^)+ij ■ 



1 



N N 

M fQ\ 



-a 



N 



a[^oo,^s<-'^ ) >0 



whence, via Lemma l3.ll . ipoo — 0. By Galerkin orthogonality for (13. 5p . a(tp — ip(p(n)i'4'ip{n)) = 0. 

II 1 1 2 II 1 1 2 

That is, ||V'<?i(n)||„ = a(V'> V'0(n))- Hence, lim„_j.oo ||V'<?i(n) ||„ = lim„_>oo a(V', '!/'</.(«)) = a{ip,tpoc) = 0. 
As the full sequence of norms (ll'/'n|la)„>o monotonic decreasing, the full sequence (V'n)„>o 
converges strongly to in H(D; M). □ 
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4.3. Rate of convergence. The theory of nonUnear approximation provides us with some esti- 
mates on the rate of convergence of Algorithm Q and Algorithm [i. Following |UT96) we introduce 
the space 

Ai := y :4°(M), (4.8) 

M>0 

where 

AliM) := e H(D; M): ^ - V cuwu, Wk e ® H(A; M,), \\wk\\, = 1, 



\K\<oo and ^|cfc|<Af|, (4.9) 



together with the norm 



:= inf {m > 0: ^ e ^?(M)} • (4.10) 

The importance of this space becomes apparent in the light of the following two theorems. 

Theorem 4.6 (Theorem 3.6 of |DT96) ). // the solution ip of (12.41) is a member of Ai, then the 
n-th error ipn of the Pure Greedy Algorithm (Algorithm]^ satisfies 

Theorem 4.7 (Theorem 3.7 of |DT96| ). // the solution ip of (12. 4p is a member of Ai, then the 
n-th error Ipn of the Orthogonal Greedy Algorithm (Algorithmic) satisfies 

Remark 6. 

(1) Pure Greedy Algorithm-based approximations such as Algorithm [l| have been proved to 
obey the slightly improved rate (see |Tem08[ Remark 2.3.11] and references therein) 

||V'„L<4||^IU^ n-"/62. 

(2) In |CEL11I Theorem 4.1] it is shown that the convergence of the Orthogonal Greedy 
Algorithm is exponentially fast if the factor spaces and the full ansatz space (in our setting 
the H(Dj;Mj) and H(D; M), respectively) are finite-dimensional. 

We note that Ai will remain the same space if in its definition — in (|4.9p , in particular — we re- 
place the energy norm || -H^ with the standard norm of H(D; M), as these two norms are equivalent. 
Then, (p S H(D; M) will be a member of Ai if, and only if, there exists an M* > such that, for 
all e > 0, there is a Xe & H(D; M) that satisfies 

lk-Xe|lH(D;M)<£. Xe=E4'^^i'^ |A(^)|<(X), 5]|ci^)|<M*; 

fceA(') fceA(s) 

and, for k € A^^), H^"^ ||h(D;M) = 1 and w^^^ e 0tiH(A;M,). 

By virtue of the isometric isomorphism described in (j2.3p . the above relations imply that 

||M-V-M-ixe||Hi(D) M-^Xs- E 4''m-^4^\ 

fceA(=) 

and, for k e Af'^), || M-iw^"' ||hi,(d) = 1 and M-^w[''' G (8>»^i Hi/. (A), the last relation being 
a consequence of the tensor-product structure of the Maxwellian M. Thus we have shown that 
M~^ip S H^(D) can be approximated to within any positive tolerance e in the norm of Hj|^(D) by 
finite linear combinations of normalized members of (S^igfAr] H^/. (-D^) with the coefficients of the 
linear combinations having their absolute sum bounded by M*. In other words, the membership 
oi If E Ai implies the membership of M~^(p in the H^(D)-based analogue of Ai, namely, 

B, := y Wm), (4.11) 

M>0 
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where 

BUM) ■■^{'pe H^D): </7 = ^ CkWk, Wk e ^ i^iu(D^), ||wfe||H^(D) = 1, 

I fcSA 

|A|<oo and ^ |cfe| < A/ L (4.12) 

feeA J 

and 

ll^llg^ :=inf{M>0: ^eSpT)}. (4.13) 

In a completely analogous way, the membership of M^^cp in Bi implies the membership of ip in 
Ai- We then have the relations 

A.^MBu IHU, - ||M-i-||g^, (4.14) 

where the last equality follows from the fact that the coefficients of the approximations to p are 
the same as the coefficients of the corresponding approximations to M~^p. 

As the definition of Ai given in (j4.8p is fairly abstract, it is of interest to have conditions in 
terms of regularity that guarantee membership in Ai analogous to the conditions provided in 
[LBLM091 Remark 4] for the separated representation strategy applied to the Laplacian defined 
on a tensor product of one-dimensional domains. This is the theme of the next section. Because of 
(|4.14p . we can pose the problem in terms of membership in the Hj|^(D)-based Bi instead with no 
loss of generality and substantial gain in succinctness; thus we shall henceforth phrase our results 
in terms of Bi rather than Ai- 



5. Characterization of a subspace of rapidly converging solutions 

5.1. Eigenvalues. We need the following two abstract lemmas, which state standard results (es- 
sentially, the Hilbert-Schmidt theorem and some of its corollaries). As we could not find these 
results in the literature in the precise form stated here, we provide brief proofs of them. 

Lemma 5.1. Let H and V be separable infinite- dimensional Hilbert spaces, with V <s H and 
V = H in the norm of H. Let a: V y, V ^ ^ he a nonzero, symmetric, bounded and elliptic 
bilinear form. Then, there exist sequences of real numbers (A„)„gjij and unit H -norm members of 
y i^n)n^fi, which solve the following problem: Find A G K and e £ iJ \ {0} such that 

a{e,v) ^ X{e,v)H ^ v e V. (5.1) 

The A„, which can be assumed to be in increasing order with respect to n, are positive, bounded 
from below away from 0, and lim„_j.oo A„ = oo. 

Additionally, the e„ form an H -orthonormal system whose H -closed span is H and the resettling 
e-nl \f\i gives rise to an a-orthonormal system whose a-closed span is V . 

Proof. This proof is an adaptation of the proof of Theorem IX. 31 in |Bre83| . The Lax-Milgram 
lemma implies the existence of an operator T : H ^ V where, given h £ H, T{h) is defined as the 
unique solution in V to the variational problem 

a{f{h),v) ^ {h,v)H VweF. (5.2) 

It also follows, via the elliptic stability estimate of the Lax-Milgram lemma and the continuity of 
the embedding V ^ H , that T is bounded. Let i: V H denote the embedding operator that 
maps V into H, i.e., v £ V ^ i{v) — v £ H . Then, T := j o T is a bounded operator defined on 
H with values in H; as i: V ^ H is a, compact linear operator, it follows that T: H ^ H is a, 
compact linear operator. Further, for all {h, h') E H x H, 

{T{h), h')H = {T{h), h')H = {h\ f{h))H = a{f{h'),f{h)) 

= a{f{h),f{h')) = {h,f{h'))H - {h,T{h'))H, 
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whence T is self-adjoint. Thus, thanks to Theorem VI. 11 in |Bre83) . there exists an i/-orthonormal 
system ien)nyi of eigenvectors of T such that 

oo oo 

/i = ^(/i,e„)_f/e„ and WHh = ^{h,en)H yheH. (5.3) 

n—l n—1 

As, for aU h e H, {T{h),h)H = a{f{h),f{h)) and a is F-eUiptic, all the eigenvalues of T 
are nonnegative. Also, as T is bounded, the set of its eigenvalues is also bounded. Now, by 
Theorem VI. 8 in |Bre83j . the set of nonzero eigenvalues of T is either empty, or finite, or countable 
with as its only accumulation point. However, on account of (|5.3p . the latter alternative is then 
the one that holds. 

If were an eigenvalue of T, there would exist e E H \ {0} such that T(e) = 0; i.e., e G Ker(r). 
However, from (|5.2p we then have that e £ V^" . As H = V ® V^" in the norm of H and V is 
dense in H , V^" = {0}, which contradicts e 7^ 0. Therefore, is not an eigenvalue of T. 

From the above, we can take the eigenvectors e„ of (j5.3[) as associated to positive eigenvalues 
fin bounded from above, arranged in decreasing order (/in+i < for ^ 1) with lim„^oo = 0- 
A consequence of the absence of from the spectrum of T is that all the eigenvectors of T have 
to be members of the smaller space V. 

Assuming that /i 7^ and e GV \ {0}, T(e) = fie if, and only if, a{e,w) ~ fi^^(e,w)H for all 
w G V. Then, all the eigenvalues of the eigenvalue problem (j5.ip are reciprocals of eigenvalues of 
T with the possible exception of 0. However, from the F-ellipticity of a, cannot be an eigenvalue 
of the problem (|5.ip . On defining A„ := fi^^ and setting the e„ to be the same as in (|5.3p we 
obtain the desired existence and distribution statements about of the eigenvalues of (|5.ip . 

We observe from a(e„,em) = A„(e„,em)fl', n > 1, that {^n/ V\i} n>i a-orthonormal 
system in V. Let us denote the a-closure of its span by V. Then, v e if, and only if, 

a{v,en) = for all n > 1. As each e„ is an eigenfunction of the problem (j5.ip associated 
to a nonzero eigenvalue, it follows from (j5.3p that v = and therefore V^^° — {0}. Thus, 
V ® V-^" = V. This, together with (lO)) itself, completes the proof. □ 



Lemma 5.2. Let the spaces H, V and the bilinear form a be as in the statement of Lemma \5.1 
and let (A„,e„) G K>o x V be the eigenpairs of (jS.ip obtained there. Then, 

00 00 
= ^(/i, e„)//e„ and = ^(^, e„)|^ \/heH, (5.4) 

n— 1 n—l 



2 



Furthe 



^ = H"("'77!=) 77^ and ||w||^ = (u, ) Vw e (5.5) 



h^H and A„(fa, e„)ff < 00 ^ e y. (5.6) 



Proof. The expression (j5.4p is just a restatement of (j5.3p in the proof of Lemma l5.lL but unlike 
there, here we emphasize that the e„ belong to V . The expression (j5.5p comes from the identity 
between V and the a-closed span of the a-orthonormal set (e„/v'A7)„>2 — P^^t of the statement 
of Lemma via, for example. Theorem VI. 9 of [Bre83| . 

As (A„,e„) is an eigenpair of (|5.ip . a{v,en/^/\[) = V^{v,en)H for all v € V; this and the 
second expression of (|5.5p give the right-to- left implication in (|5.6p . Let us now consider a.n h E H 
that satisfies the left-hand side of (|5.6p . As the e„ are members of V, the partial sums 

fc 

hk ■= ^^{h, en)Hen 
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also belong to V. The a-orthonormality of the e„/vXn leads to the equality, for I < k < I, 



I 

2 



\\hi-hk\\l= ^ A„(/i, e„)^. 

n=k+l 

As the real series J2^=i ^n{h, e„)|^ is assumed to converge, the above expression tends to as fc 
and I tend to oo. Hence, (^fc)fe>i is a Cauchy sequence in V and thus converges to some h V. 
As V is continuously embedded in H (part of being compactly embedded), the limit h has to be 
the same limit the hk have in H. That is, h = h G V. This completes the proof of (|5.6p . □ 



The hypotheses of Lemma |5.1| and Lemma |5.2| are satisfied by the eigenvalue problems 

{e^'\v)Hl,SD.) - AW(eW,^)L.^^(^,) e H1,^(A), (5.7) 

(for i E [N] here and in what follows), and 

(e,¥')H^(D) = A(e,^)L2(D) V(y3eH^(D), (5.8) 

whence their solutions do have the distribution, orthogonality and spanning properties stated in 
that lemma (the hypothesis V = H , which is not discussed elsewhere, follows from the density of 
infinitely differentiable and compactly supported functions in any weighted L^ space). In particu- 
lar, they have sequences of solutions (eigenpairs) ({Xn\en^)j and ((A„, e„))^gpj, respectively, 
with 

oo 

VeLluiD,) and ^ (^, e«)2,^ (^^^ < oo ^ ^ G H],,, (A), (5.9) 

n=l 

and 

oo 

^eL2,(D) and ^ A„((/7, e„)22 < oo G H1,(D). (5.10) 

Next, we exploit the special tensor-product structure of the full Maxwellian M to characterize the 
eigenpairs of its associated eigenvalue problem (15. 8p in terms of the eigenpairs of the eigenvalue 
problem (j5.7p associated to the partial Maxwellians Mi. 

Lemma 5.3. The net ((An, e„))^^j^„^ „„)gNN is a full system of solutions of the eigenvalue 
problem (|5.8p . where 

A„ := 1 + ^(AW - 1) and e„ := e^. (5.11) 
Proof. Given r = t(*) G ^ie[N] '^o'l^i), we have that 

N N 
AT N 

Since the span of (^^^ H(Z?i; Mi) is dense in H(D; M) (as is readily seen from Corollarv l4.4l and 



(12. Sp ). the equality of the first and the last expression in the chain of equalities above is valid for 
all r G H(D; M). Hence, (A„,en) is an eigenpair of (|5.8p . Further, we deduce from the chain of 
equalities above that is orthogonal to Cm in both L^(D) and H^(D) if n 7^ m. 

From ([53)) in Lemma [siil. for i G [N], span (e^*^) = R\jXDi). Hence, by (03) and 

V / n>l ^ 



(g) span (e^n^) C span (e„)„gpjjv = Hj^(D). 

i=l ^ ) n>\ 
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Thus, (e„)^gpjN forms an orthogonal system that spans H^(D). Therefore, by Theorem VI. 9 
of |Bre83| . all the eigenpairs of the (full) Maxwellian eigenvalue problem (jS.Sp have the form 
(An,e„) as given in (jS.lip (modulo linear combinations of eigenfunctions belonging to the same 
eigenspace). □ 



It follows from Lemma l5.3l that the eigenvalues and eigenfunctions of (j5.8|) are more naturally 
indexed by than by N; in what follows, we shall refrain from indexing contra natura. 



5.2. Characterization via summability of Fourier coefficients. As b y Le mma l5.3l the se- 
quence ((A„, en))„gpjiv is a full system of eigenpairs of (|5.8p . (|5.5p in Lemma |5.2| ensures that, for 
all re Hi, (D), 

= in Hj^(D). 




Hence, given the tensor-product structure of the and the unit H^(D)-norm of the e„/%/ we 
can guarantee that r e Si (cf. (|4.1ip ) if 



< OO. 



In turn, this holds if 



A. 



A:— ^ — < OO and B := ^ o-„(r, eri)L2 (□) < oo, (5.12) 

where ((T„)^gf^jv is a sequence of positive real numbers that are to be chosen below. We note 
that the requirement of B being finite can be seen — for (t„ — A„, for example, this is certainly 
the follows from (|5.10p — as a regularity requirement on r. Thus, there is a trade-off in 

(I5.12[) between the requirement that the (T„ grow fast enough to ensure the finiteness of A and the 
desirability of the an growing slow enough to avoid demanding more regularity than necessary of 
the functions r for which B is finite. 

As a first step in formalizing the above we consider, given a net S = (cn)„gNiv with entries in 
M>o, the space of all those L^(D) functions for which the term B, as defined in (|5.12p . is finite: 

H^(D) := L e L^(D) : ^ a„ e„)^. < oo I . (5.13a) 



Tie: 



IN 



We equip Hy(D) with the norm 

\ 1/2 



II'^IIhS(d) ■= I ^ f^"(¥',en)L2(D) J • (5.13b) 

\neN« / 

It is readily seen that, if there exists a cr > with > ct for all n G N^, then H^(D) is a 
separable Hilbert space that is continuously embedded in L^(D). Further, if there exists a cr' > 
such that On_> <^'^n for all n e N^, then H^(D) is continuously embedded in H^(D) and, thanks 
to Lemma [2!^ it is compactly embedded in L^(D). 

At this stage we could just choose S to be, e.g., cr„ — Xn \\n\\2 for some a > N and an 
application of a multiple series version of the integral test for convergence (see, for example, 
[GLlOi Proposition 7.57]) would render the sum A in (|5.12p finite. However, the resulting space 
H^(D) would then still have quite an abstract description. What we therefore wish to do instead 
is to choose each ct„ as a suitable polynomial function of the {X^^\ . . . , A^^-'). Then, under certain 
reasonable conditions, which we will make explicit below, we shall be able to characterize the 
resulting space in terms of regularity properties. One of these conditions has to do with the fact 
that we can only know that A of (I5.12p is finite, with certain polynomial function of the 

An, if we have some information about the asymptotic behavior of the An. Consequently, we adopt 
the following hypothesis. 
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Hypothesis C. For each i e [N] there exist positive real nurabers cf' and and n^*) e N such 
that (with d signifying the common dimension of the single-spring configuration domains Di) 

where An'' is the n-th member of the (ordered, with repetitions according to multiplicity) sequence 
of eigenvalues of (|5.7p . 

Remark 7. Hypothesis basically consists of assuming that the eigenvalues of the problem (|5.7p 
behave like the eigenvalues of a regular elliptic operator, such as the Poisson operator. If the partial 
Maxwellian Mi com es fr om either the FENE model ([12]) or the CPAIL model ((O]) Hypothesis 
holds; see Corollary C.2 in Appendix for a proof (see also Remark 11). 

Theorem 5.4. Let T^™) = (^t4"'') be defined by 

N 

ri"^-n(^nO" VneN^. (5.14) 

i=l 

Then, H^'"' (D) C 6i i/ m > f + 1. 

Proof. According to the previous discussion, the stated inclusion will hold once we have shown 
that the infinite sum over n G of Xn/rn^'' converges; i.e., that A in (|5.12p is finite. To prove 
this, we start by noting that, modulo a decrease of c\ and an increase of , we can take n^'-' — 1 
in Hypothesis as a consequence of all the Xn^ being positive; we do so from now on. This, 
together with (|5.1ip and Hypothesis [c|, yields that 

for all n e and some C > that depends on the c^^\ the C2*\ N and d only. Clearly, it will 
be enough to show that the right-most expression in (I5.15P results in a convergent series. Now, 

( N 00 ^ 



v-^Af 2/d N 2/d-2m/d N 

sr^ l^i=\ "i - 'S^ 'S^ 3 ^ 

2-^ 2-^ T-iN 2m/d 2-^ 



00 

n E 



J2 



=1 



where the constraint on m ensures that all the resulting one-dimensional sums are finite 



□ 

For later reference we introduce another family of weights that also produces subspaces of B\. 
Theorem 5.5. Let T^"') = { v~f^\ be defined by 

/ N \ ™ 

^i"^:=(E^«'j VneN^. (5.16) 

Then, HJ5,*'"'(D) <Z Bi if m > 1 + \Nd. 

Proof. Using Hypothesis [c| and the already mentioned multiple series version of the integral 
test for convergence it can be shown that the result hinges on the finiteness of the integral 

/[I 00)" (X^iLi ^^i^'') dx. Thanks to the equivalence of the 2/(i-quasinorm to the 2-norm in 

^>0- 11^112 

turn, implied by the finiteness of the integral 



Euclidean space and since [l,oo)^ C {a; G K>q: ||a;||2 > 1}, the finiteness of that integral is, in 



x\ 



9 / ^ 
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where Cn is the (A^— l)-dnTLensional volume of the surface {x £ R>q : ||a;||, = 1}. As it is assumed 
that m > 1 + ^Nd, the last of the above integrals is finite and the proof is completed. 

□ 

The definition of H^""' (D) given by (|5.14p is less abstract than the definition of Bi (given in 
(|4.1ip '). However, we can describe subspaces of the former space in even less abstract terms by 
showing that certain regularity conditions translate into summability conditions expressed in terms 
of Fourier coefficients, such as those that define H^'"' (D) (cf. (I5.13ap ). In order to understand the 
appropriate regularity requirements for this purpose, we need to study the regularity properties 
of certain degenerate elliptic operators in Maxwellian-weighted Sobolev spaces. 

5.3. Characterization via membership in mixed-order weighted Sobolev spaces. We 

start by adopting two further hypotheses. 

Hypothesis D. For i £ [N] the spring potential Ui is monotonic increasing and convex. 

Hypothesis E. For i £ [N] there exists a distance ji £ (0, \/bi), an exponent > 1 and a 
function hi £ ([0,7i]) that is positive on [0,7^], such that 

Mi{p)^h,{Mp))Mpr 

for all p £ Di such that i)i(p) £ (0,7i), where 0^ is the distance-to-the-boundary function in Di. 

Remark 8. Hypothesis [3 can be regarded as a strengthening of Hypothesis 0- It is easy to check 
that springs obeyingthe FENE model (jl.2p or the CPAIL model (|1.3p comply with it. 

With Hypothesis |E| we are restricting ourselves, essentially, to power weights. The compliance 
of the FENE and the CPAIL models with it is also easy to check if their parameter bi is greater 
than 2 in the FENE case and greater than 3 in the CPAIL case. 

Lemma 5.6. For i £ [N], 

(a) the space C^{Di) is dense in }i\j,{Di); 

(b) the space C°°(A) is dense in H^}^(A), for m € N. 

Proof. In Proposition 9.10 (resp. Theorem 7.2) of |Kuf85| the result (a) (resp. (b)) is stated for 
weights that are powers greater than 1 (resp. greater or equal than 0) of the distance-to-the- 
boundary function; the bilateral boundedness of the function hi by positive constants, implied by 
Hypothesis 0, extends the statement to our case. □ 

The additional requirements on the potentials Ui, i £ [N], and the preceding lemma allow us 
to prove a first elliptic regularity result. 

Lemma 5.7. Let i £ [N] and suppose that g £ \^\^,{Di); then, there exists a constant Ci > 0, 
independent of g, such that the solution z £ Jl\^j\Di) of 

{z,^)bI,Ad.) {9,^)Ll,^iD,) £ H],,j^(A) (5.17) 

obeys the regularity estimate 

1 



Mi 



Proof. By Hypothesis and Hypothesis the function : — > M defined by Vi :— |A(^ M^) 
is convex and tends to +oo as its argument approaches the boundary of Di from within. Then, it 
follows from Theorem 3.4 of |DPL04| and the density of Cg°(Z3i) in H^.(Di) given in part (a) of 
Lemma [53 that there exists a unique solution z in {u £ H^j. (D^) : VVi • Vm G L\j.{Di)} to 

iAz + VF, • V5 = ig, (5.18) 
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considered as an equation in L^j.(Z3i), and it obeys the estimates 



(5.19a) 
(5.19b) 
(5.19c) 

The regularity of Mi and z admits the use of the Leibniz formula for t he p roduct of a regular 
distribution and a continuously differentiable function provided in Lemma IB.II in Appendix |B|. We 
can then write M^Az — 2Mi\JVi ■ V z = div(MiVz) (for this we have used that Ali is proportional 
to exp(-2T/,) (cf. dUS])). Plugging this into ([ETS|) gives 





\\9\ 












= ^\\9\ki^XD 








= 2|l3llLi,^(D.) 



Mi 



div(M,V5] 



< 



Li.. {D. 



(5.20) 



Multiplying (|5.18p by 2Mi and using the Leibniz formula for the product of a regular distribution 
and a continuously differentiable function again, we find that 



zipMi 



Di 



Vz • VifMi 



Di 



9V> 



for all (fi e Cg°(Z)i). It follows from the density of C'^{Di) in {Di) and the uniqueness of the 
solution z of (j5.17l) that z = z and hence ()5.19p and (j5.20l) give the desired result. □ 

In order to obtain an iterated elliptic regularity result, we need the technical lemma that follows. 
Lemma 5.8 (Hardy inequalities). Let H > 0. Then, there exists C'h > such that 



v/eLi((o,ff)). 



(5.21) 



If a> I, then there exists CH,a such that 



H 



H 



y"-'f{yfdy<CH,a y'' [f {yf + f iyf] dy V/ e 4^4(0, i/)). (5.22) 
Jo 

Proof. The inequality (|5.2ip follows from the standard Hardy inequality (the H — oo case); see, 
for example, |DiB02[ Proposition VIII. 18.1]. Alternatively, see jOK90i Theorem 1.14] for a very 
general form, which encompasses (I5.2ip . 

To prove (|5.22p we will use a procedure inspired by the proof of Theorem 8.2 of [Kuf85j . The 
first ingredient is the inequality 



r y'-'fiyrdyKC^ [": 

Jo Jo 



y'^nyfdy 



valid for aU / in G^{[{),H]) such that f{H) = (see, e.g., [OK90[ Example 6.8.ii]). Let now ip^ 
and Lpi form a smooth partition of unity subordinate to the covering H = (0, 2H/3) U {H/3, H). 
Then, given any / G C^([0,i/]), let /o ■— fof and /i ipif. Using the above inequality, the 
vahdity of ((O^ for C^{[0,H]) functions follows from 



lL^.,<,_2((0.ff)) - ll/ollL^^^„_,((o,2ff/3)) + II/: 



< Cl^^ II/oIIl^,j„((0,2H/3)) + ||(") ^M\lI^^{(H/3M)) 

< Cl^^ \\(pof' + <^o/|Il2^„((0,2H/3)) + II/iIIl2j„((H/3,H)) 

< (^2 II/'IIl2 ((0,2H/3)) + II/IIl^ ((0,2H/3)) + II/IIl^ ((ff/3,H)) 



+ \\r\\ 



Lf.)„((0,H)) 
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The validity of the inequaUty for all / G Hj ((0, i?)) is then a consequence of the density of 
C^([0, H]) functions in ((0, H)), the completeness of .|„_2 ((0, H)) and the continuity of the 

2 



injection of that latter space into L( \a ((0, H)). □ 



We shall now iterate Lemma |5. 71 : extra regularity for g implies extra regularity for z. 

Lemma 5.9. Let i g [N] and g £ H^j. {Di). Then, there exists a constant Ci > 0, independent of 
g, such that the solution z G H\,j,{Di) of 

{z,^)Hl^iD,) = {9,^)lI,AD0 V(y5 e H]v,^^(A) (5.23) 
obeys the regularity estimate 

ll^lln^j^Di) - llffllHljJDi) • 

Proof. The core of this proof is based on Lemmas 3.1 and 3.3 of |Fre87j . As their adaptation to 
our geometry is nontrivial, we give a detailed argument. Note that in this proof we shall omit the 
spring index i in order to avoid cluttering the notation. 

Part 1 : We start by describing a change of coordinates and how (j5.17p transforms under it. 

Given p e M.'^, let p' denote {pi, . . . ,Pd-i) G M''^^. Let ( be some constant in (0, 1) and let us 
define, for e £ (0, (], the sets 11^ := x (0, £7) and [4 :— S{Ue), where 7 is the distance (with its 
spring index omitted) mentioned in Hypothesis [2 and 

p, ,^ ie{--/2,^/2) ifd = 2, 

" ' |e(-l,l) X £(-V2, V2) ifd = 3 

and 5*: C/^ — >■ C/^ is defined by the formula 

f(V^-P2) (cos(pi), sin(pi)) ifd = 2, 

\(N/6-p3)(yT^cos(p2), yr^sin(p2), Pi) ifd-3, 

Note that if < £1 < £2 < C then [/^^ C 11^2 C D and C/^^ C C/e^. Having its domain and 
defining formula carefully crafted for the purpose, the transformation S turns out to be invertible, 
orientation-preserving and C°°(c7^)-regular. All of this is easy to see if one takes into account that 
S" is a variant of the polar (resp. spherical) to Cartesian coordinate transformation if d = 2 (resp. 
d = 3) with the radial variable being measured from the boundary of D and increasing towards its 
center. We denote the inverse of S by T; it, too, has uniformly bounded derivatives of all orders. 

If / is a function with domain ?7<^ we will write / := foS. Then, ip e C°°{U(;) <J=^ € C°°{uc)- 
If m is a positive integer, part (b) of Lemma l5.6l states that C°°{D) is dense in H^,}(Z?); as, for any 
£ G (0, C], Us is regular enough (a Lipschitz domain), C°°(?7e) is exactly the set of restrictions t o Uf 



of C°°{D) functions, whence C°°{U^) is dense in H^(J7e) as well. We also have from Lemma lB.3 



in AppendixlBl, that / e H™ (t/^) / G B^^^iUe) and that 

ci(™)II/IIh;.(c/.) < WfKiUU.) ^ C2(m)||/||H™(^^) V/ G R^iU,), (5.24) 

where the positive constants ci and C2 depend on m but can be chosen to be independent of £. 

As C°°{U(^) is mapped by composition with S to C°°(i/j) bijectively, it follows that C°°(c/<;) 
is dense in H^~^([/^). The rules of calculus and the density of C°°(J7,j) and C°°{U() functions in 
■'^m(^c) ^^'^ "'^A/^ I'espectively, give the equalities 

/ uvM ^ / uvMa and / Vu • VwAf = / VuAVv^M, (5.25a) 
Ju^ Jiff 



where we have used the shorthand notations 

a = det(VS') and A = {V S)-^{V S)''^ det{\/ S). (5.25b) 

The first equality in (j5.25a[) is valid for u and v in h\j{U(^) and the second for u and v in H]^,^(J7^). 
Direct calculations give 

A = A^, Ak,d^Ad,k=0 Vfce[d~l], Ad,d^a, and dka^O Vfce[d~l], (5.26) 
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Figure 5.1. Illustration of the construction used in the proof of Lemma Is.Ol . 
From left to right: Contour plot of M on D with [/^ and enclosed by the 
thick continuous and dashed lines, respectively; contour plot of M on Uq with Ug 
enclosed by the thick dashed line; contour plot of an admissible uj on Uq with 11^ 
enclosed by the thick dashed line 



which we will exploit later. The need for the last equality in (|5.26p is the rationale behind 
taking of the sine of the polar angle instead of the polar angle itself as the first argument of the 
transformation S in the case of d = 3. Additionally, by construction, M is a function of the radial 
variable pd only — namely, for all p S f/^, M{p) = h{pd)Pd, where h and a, with the spring index 
omitted, are those of Hypothesis 

Let us fix e g (0, C)- For localization purposes we pick a C°°(C/^) function oj with range [0, 1], 
identically 1 in U^, with support bounded away from OUq \ dD and such that ddio{p) = if p is 
within a finite distance 7' > of dUi^ n T{dD) (such a function is readily constructed as w = (D o T 
where w(p) = s{p')t{pd) and s and t are suitable mollified step functions). See Fig. IS.ll for a 
depiction of the construction so far. 

Now, as every member of Cg°(f/(;) can be put in the form foS, where (p G Cg°(C/^) C H]y,j(J7^), 
(|5.23p and the equalities in (|5.25p imply that z obeys the distributional equation 



div (VzAAl'j + zaM = gMa in [7^. 



(5.27) 



Part 2 : In this part we show that the relevant derivatives of z in directions tangential to 
the radial (i.e., the d-th) coordinate possess additional regularity. The argument is a nontrivial 
adaptation of Lemma 3.1 of |Fre87j . 

Let fee [d — 1] . Then, using the Leibniz formula (cf . Lemma iB.lh and simple consequences of 



(I5.26|) and the fact that M{p) depends on pd only, the distributional equation (|5.27|) conduces to 

- div (y {ujdkz) AIvI^ + ujdkz aM 

= Cjdk9aM + LbWz : duAM + CjV z ■ d\v{dkA)M + CoV M ■ {VzdkA) 

- TsJCj ■ (Vidkz) AM) - dkSdiv{VCjA)M - dkzVM ■ (Vw A). (5.28) 

We want to show that all the terms resulting above are (the linear combination of) members 
of the space aML^^([/|^) = ML^-^([/^). Of the resulting seven terms above, the first three, the 
fifth and the sixth pose no problem, thanks to the regularity of g and Lemma |5.7| . The fourth 
vanishes after making full use of the equalities in (|5.26p and the sole dependence of M on the 
radial variable — this is what the fourth equation in (j5.26p is truly for. The membership of the 
seventh term in ML'^- (Uq) stems from observing that 



dkzVM ■ {VCjA)/M 



dkZ ddM ddCo a/M 
sup 



< sup (ddM ddUJa/M] |15a;z|1l2 (fj. < 00. 

P^x(7',C7) ^ ^ " ' 



Let /, given some function / defined on U(;, denote f o T = f o S ^. Also, let f/^k) denote the 
ratio of the right-hand side of (I5.28P and Ma. Then, (j5.28p and the identities in (|5.25p give 



div(MVa)9/jz) + CodkzM — f(^k) 



(5.29) 
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in Uq, with /(j,) G L^,j([/^). As the support of u is bounded away from dUi^ \ T{dD), we can 
extend uidkZ and /(fc) to the whole of D by zero while still satisfying (j5.29p . Then, Lemma IsTtI 
ensures that the extension of ujdkZ belongs to H\j{D). It follows that QdkZ belongs to H^^(J7(;) 
and, consequently, 9^2 G H^-^([/e). 

This procedure can be iterated. Within the identity (|5.28p particularizes to 

-div (v(9fcz) AM^ + dkzaM = dkgaM + VV5 : dkAM + ■ div{dkA)M. 

Let g(^k) ■= dkg + VVz : /a + Vz • div(9fcA)/a and let us redefine ui so that the role of is 
now taken up by and the role of the latter is taken up by Us, where S is some fixed number in 
(0,£). Thus, we can obtain an analogue of (|5.28p for ujdi^kZ, where l,k G [d— 1]: 

- div {\7{ujdi^kz) AM) + LddLkzaM 

= Cjdig(k) aM + uWdkZ : diAM + CjV dkZ ■ dhj{diA)M + wVAf • (V(5fcz) diA) 

- 2VCj ■ (V(ai,fez) AM) - di^kS div{y Cj A)M - di^kSV M ■ (VtiA). 

Analogously to the study of the first-order tangential derivatives we need all seven terms on the 
right-hand side of the above equation to belong to ML^(C/e) now. As, at this stage, we know that 
dkZ £ H^-^(C/e), the second, the third, the fifth and the sixth term above pose no difficulties. The 
fourth term and the seventh term can be dealt with just as their counterparts in (|5.28p . When it 
comes to the first term, it is enough to show that dig(k) G L^~^([/e). Now, 

dWik) = diyg + VV~z : diidkA /a) +\7 z ■ diidividkA)/a) + diW 2 : {dkA /a) + di\7 S ■ div {dk A) /a. 

The first three terms above are clearly in L^(C/e) — the first because of our hypotheses on g; so is 
the fifth, for the second derivatives of z have the desired integrability. Finally, diS/S/z G L^~^(C/e) 
because diz G H^-^([/e), as shown above. Proceeding with the argument one finds, after localization, 
that dk,iz G H^-^(t/5). We mention in passing that by closely following the arguments above the 
linear operators g G ll\.j{D) i-> diz G H^-^(f/e) and g G H^j(D) dkjz G ll^^~^{tJs) can be seen to 
be continuous; i.e., bounded. 

Part 3 : In this part we show the additional regularity of some derivatives of z that involve the 
radial direction. 

Expanding and rearranging the distributional equation (j5.27p . noting the sole dependence of 
M on the last component of its argument and the properties of A given by (|5.26p we get 



-^ddiddzaM) = g- z + - ^^{dj^kZ Aj,k + djzdkAj^k) =: / (5.30) 
Ma a , — ■ 

k=l j = l 

in Us- From the previous part of the proof and our assumptions on g we have that / G Y(^^^{Us)- 
Multiplying (j5.30p by Ma and integrating with respect to the c?-th variable we obtain 

rPd 

{ddzaM)[p\pd]^ \iuY{dd~zaM)[p',s]^ {faM)[p\s\ds (5.31) 

for almost every p' in P'g. We note in passing that in this part of the proof we reserve square 
brackets for arguments of functions. Our first task is to show that the limit on the left-hand side 
of (|5.3ip vanishes. To this end, we first observe that, for pd, s G (0, ^7), 



q/2 
Pd 



— da. 



whence 



Pd 



/2 



\dd~z[p',Pd]\<s"/^\ddS[p\s] 



/Pd ^ rPd 
-a"/^-'dd~z[p',a]da + J a''/^dd,dS[p' ,<j]da 
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Furthermore, 



Pd \ddS[p',Pd]\ 

<3s"|adz[p',s]| 



<3s"|9dz[p',s]r 



4 
4 

„2 



/■pd ^ /"Pd ^ 

/ a"/2-ia^|[p>]da +3 / a''/^dd,dS[p',a]da 

J s J s 

fPd rPd 

\pd-s\ a''-^\ddS[p',a]f da + 3\pd~s\ a" a] |' da 

J S J S 



<3s"\dd~z[p',s]\^ + ^d^ r a"-^\dd~z[p',a]f da + 36^ H ^" \dd.d~z[p' , da. 
^ Jo Jo 

Integrating this chain of inequahties with respect to s from to ^7 and applying the Hardy 
inequahty (|5.22p stated in Lemma IS.Sl we obtain 



^7 



S-/p2 \ddS[p',Pd]\ 

nS'y 



< 3 



s'^\dd~z[p\s]f ds+^{Sjf 



r a^~^ \dd~z[p\af da + 3{6-if /'^a" \ddA~z[p\'yf da 
Jo Jo 



< 



3 r s'^\ddS[p',s]f ds+ 5^^^(57)2 r a'^\dd~z[p',a]f da 
Jo 4 Jo 

+ 3(^^)2 (^^!^ + 1^ \dd^,S[p',af da. 



Dividing by Sj, integrating with respect to p' in P^, using the bilateral boundedness of h and a 
by positive constants, and consolidating the constants, we get the trace-inequality-like bound 



p'^\dd~z[p',Pd]\ dp' <Ci\\ddS\\li^ 



(Us) 



(5.32) 



Thus, 



Pi 



iddSaM)[p',pd] dp' < C2h[pdV^'p'^ 



1/2 a/2 



1/2 



Pi 



UddSfMaMdp' 



as Pd — > 0+, 



which implies the vanishing of the limit in (|5.31[) . 
Let us define w: t/^- — > R by 



dd(Ma)[p] 
{Ma)[p] 



ddz[p] = 



iiha)[pd]p2y r 



iha)[pd]^pj- 7o 



{faM)[p',s] ds, 



(5.33) 



where we have taken the liberty of treating a as an univariate function, which it is in the algebraic 
sense. The equality is valid for almost every p € IJs. Note that is a member of L'^^j{Us) because 
VM • Vz/M G L\.i{Us); this, in turn, is a consequence of Lemma 15.71 . We intend to show that 
w G H^([/5). Let (/n)„gN be a sequence of C°°{us) functions converging to / in R'^^^{Us) (its 
existence having been discussed in Part 1) and let 

{{ha)[pd]p2)' r, 



Wn[P\ 



{ha)[pdYpl" Jo 



ifnaM)[p',s] ds 



{ha)'[pd]pd + a{ha)[pd] ( s \ " {hafn)[p' , s 



{ha)[pd] 



Pd 



Pd 



■ ds = aux[prf] / S.'^{hafn)[p' ,Pd^] d^ 
Jo 



(5.34) 



Here we have written aux[prf] in place of the first fraction in the second integral and denoted the 
function p G i-> h{pd) G M by ft- as well. The second equality comes via the change of variable 
^ = s/pd- As the function h and the determinant a have uniform C'^ and G°° regularity in Us, the 
function aux G G'^{us) and Wn is twice continuously differentiable in the d-th. direction. 
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Differentiating the last integral representation of w„ with respect to its d-th variable twice and 
then reversing the change of variable we obtain 



dd,dWn[p] = ^ 



fc=0 
2 



2\ d"^ '^aux. 



-AkJ dp- 



E 



2\ 9^ '^aux , 



-oVM dp: 



2-k 



[Pd] / d^{haU)[p',s] 



a+k 



■ ds, 



whence, as s/pd € (0, 1) if s G {0,pd), 

2 



Pd^^ \dd,dWn[p] \ <— I [Y] 

Pd Jo 



Pd 



< 



1 



Pd 



< 



Pd Jo 



k=0 



E 

\k=0 



d"^ '^aux 



dPd 
i9^~''aux 



-[Pd]a,^(/.a/„)[p',s] 



Pd 



a/2+k^ 



dpi 



— \Pd]d2{hafn)[p',s\ 



\{haUf + \dd{haUf + \ddAhafn)n [p' , s\s''l^ ds 



1/2 , 



Pd Jo 

for some C3 > independent of p = {p' ,Pd)- We square the resulting i neq uality, integrate it with 
respect to pd from to (57, use the Hardy inequality (|5.2ip in Lemma IS.Sl and note yet again the 
bilateral boundedness of h and a by positive constants to obtain 



(57 



p'^{ha)[pd] \dd,dWn\p\\^ dpd 



<Ci r (KWrOl" + \dd{hafr,)f + \ddAhafn)f)[p] s"(HM dprf, 



where C4 is still independent of p' G Pg. Integrating this with respect to p' e Pg, using the 
regularity of h and a and taking into account that (Ma)[p] = {ha)[pd]p'^ for all p E Us one gets 



\\dd,dw. 



< C's ||/n||H2 



(Us) 



This argument can be carried over to all derivatives of order less than or equal to two of w„ 
(including zeroth order derivatives of w„, meaning w„ itself). The result is 

ll'^"llH2.^(C/a) ^ ^6 II/«IIh2^((75) • 

As H^(t/5) is a Hilbert space, there exists a subsequence (^^0(n))„>]^ with a weak limit w* E 
tP^{Us). By the continuity of the injection of lP^{Us) into L'^^^{Us), w* is also the weak limit of 
the in ^'^^{Us)- 

Now, given any x e L^^(?^5), 



Us 



au-x[pd\ 



Pd Jo 



Pd 



{hax)[p\s]ds] M[p]dp 

aux[prf]^ 
Us Pd 



Pd 



\Pd 



a/2 



W2(/iax)[p',s]ds h[pd]dp 



< 



f 1 / f^"^ \ 

/ - / s"/2(/iax)[p',s]ds dpddp' 
"'P^ "'0 Pd \Jo J 

Cs [ r s''\ihax)[p',s]f dsdp' 
J Pi Jo 



< C9 Wxwi^ 
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Hence, the operation that defines w (rasp. w„) in terms of / (resp. /„) in (j5.33l) (resp. (15. 34^ ) is a 
bounded map from L'^^{Us) to itself. Therefore, hm„_j.oo fn = / in L'^^j{Us) imphes hm 
w in the same space. Thus, w and the weak Hmit w* have to be the same measurable function and 
so u; g H2^(J75). We get the bound II w||h2_ (j^^) < liminf„^oo|k^(„)llH2^((75) ^ C!6\\f'l>{n)\\Hl^(Us)- 
As (with no loss of generality) we can assume that the /„ are scaled so that their lP^{Us) norm 



is identically equal to the same norm of /, it follows that 

< Ce II/IIh^^ 



w\ 



From (jOO)) and ((OSl) . 



-dd,dz / + 



ddjMa) 
Ma 



(Us) ■ 



ddz^ f - 



w, 



whence ||9d,d^|lH2 (fj,) < (1 + Ce) \\f\\^2 ,(j ) < Cio \\g\\ii2 ,jj y We know from the previous part 
that all second derivatives of z that do not involve the d-th direction have H^-^(f75) norms bounded 

by the H^^(t75) norm of g. This and the corresponding result for dd,dZ is enough to be able to 
bound all derivatives of z of order less than or equal to four, and thus deduce that 

or, equivalently in the light of (j5.24p . that 

M^i,{u,) < MhUUs) ■ (5-35) 

Part 4 : By modifying the transformation S one can get a localized bound of the form (j5.35p for 
any origin-centered rotation of Us- It follows that (|5.35p remains valid (with some other constant 
C12) if we replace Us by the annulus/spherical shell {p £ D: \p\ > Vb — Sj}. 

Let Do be the baU B{0, Vb - S-f/2) d D. As C^{Do) c Cff (L»), we have that 

The existence of a positive infimum of M in Dq implies that z is the weak solution to a regular 
(i.e., uniformly) elliptic problem in Do with H2(Do) right-hand side. It follows, via the C^'^{Do) 
regularity of M (see, e.g., |GT01[ Theorem 8.10]), that for some C13 > 0, 

ll^llHt,(D^) ^ <^13 Mul^iDo) ^ '^13 \\9\\hI(D) > 

with D'o B{0, y/E - 3(57/4) d Dq. 

Combining this last estimate with the result in the annulus/spherical shell mentioned above 
(which in union with Dq covers D), we obtain the desired global bound. □ 



The next lemma is an almost trivial corollary of Lemma l5.9l yet it is a true iterate of Lemma l5.7 
in the sens e that the hypothesis on the right-hand side function is the thesis on the soluti on in 
Lemma 5/7. This makes it suitable for the arguments that will be used in the proof of Lemma Is. 111 . 



Lemma 5.10. Let i e [N] and suppose that g e ll\j.{Di) and that ^div(A/iV.g) G L^iX^i 
Then, there exists a constant Ci > 0, independent of gi, such that the solution z G (Di) of 

obeys the regularity estimate 



— div(Af,Vz) 



-div M.V 



— div(Af,Vz) 
IVl,, 



(DO 



<c. 



'■9I1hL(d.) + 



— div(M^Vg) 



Li.. (-DO/ 



Proof. This follows directly from Lemma Is.Ol on noting that ^div(MiVz) = 5 — z in the 
distributional sense first, and then in the sense of measurable functions. □ 
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Lemma 5.11. Let i G [N]. The following statements of equivalence hold: 
reH^^^(A) and ^div(M,VT) e L^^^ (A) 



reLi,^(A) and ^(A«)2(r, eW)^, ^^^^ < oo; (5.36) 



n=l 



and 



reH4,^(A), -;Vdiv(M,Vr)eH2,,^(A) and -3-div ( Af,V 



— dW(MiVT) 



reLi,^(A) and £(A«)4(r,, 



< cx). (5.37) 



Proof. Wc will omit the spring index when proving (|5.36p and (j5.37p . We start by denoting by 
L the operator that associates each cp G W[q^(Z?) with the distribution A/~Miv(AfV(^) (this is a 
well-defined distribution because of the regularity of ip and M; cf. Lemma Ib.iI in Appendix [sh. 
We also write L := —L + I where / is the operator that associates each distribution with itself. 
Let us define the Hilbert spaces (and associated norms) 



Rl,{D) := e H|,(Z?): L{^) e LI,{D)} 



and 



Mli,{D) ■= II^IIh1,(d) + ll^(<^)lll^„(Z3) ■ 



(5.38a) 
(5.38b) 

(5.39a) 
(5.39b) 



Because of the definition of H^^(-D), L: H^^(Z?) — > L\,j{D) is a bounded linear operator. As for 
every £ L\,j{D) the solution z £ R\,.[{D) to 



exists and, thanks to Lemma l5. 71 . is bounded in R\,j{D), L ^ : L\j{D) R\^{D) is well-defined 
and bounded. Similarly, by the definition of R\j{D) and Lemma 5. id : R\.[{D) — > L\.j{D) is a 
bounded linear operator with a bounded inverse. 

Let T g R\.j{D); i.e., r complies with the left-hand side of (j5.36p . It then follows that fr 



L{t) + t £ Ij\AD) and Parseval's identity thus yields 



oo > WfA 



ri>l ri>l 



'^«/L|,(D)' 



n>l 



where (/r, e„)L2j(i3) = (•r, e„)Hij(i3) follows by the density of G'^{D) in R\AD). 

To prove the converse, note that the eigenfunctions e„ of (15. 7p are solutions of e„ — i~^(A„e„), 
whence ||e„||g2 ^^-j < C || A„e„||L2 (-q-) = CA„. Consequently, the partial sums 



n=l 



are members of H|,j(D). Then, if fc < the L^j(_D)-orthonormality of the e„ yields that 



L(rO - i(Tfc) 



As the sum Y.n>i ^li'^^ ^n] 



(■^''S")Li„(r>)i(en) 

ri=fe+l 



is assumed to converge, the sequence (i(Tfe)) is a Cauchy 



sequence in L1j{D) and hence it converges to some /* £ Ll,j{D). The continuity of L ^ implies 
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that the sequence {Tk)f.^i converges in H\.j{D) to L ^(/*). The same sequence converges in L\j{D) 
to T. The continuity of the injection of H\j{D) into L\j{D) then imphes that r = L^^{f*) E 
lilj{D). This completes the proof of ((5736t . 

Let us suppose now that r in Hj^^(_D); i.e., r comphes with the left-hand side of (j5.37p . It 
follows that fr ■= —L{t) + t e H|,^(£>) and gr ■= —L{fr) + fr ^ L^^(£)). Parseval's identity gives 



n>l 



n>l 



n>l 



n>l 



where the second equality follows, similarly as above, by the density of Cjf (Z?) in ]i\,j{D) thanks 
to the boosted regularity of fr- The latter also allows the use of (|5.7p to obtain the third equality. 

To prove the converse we first note that each e„ is a solution of e„ = L^'^{\\en), whence 
ll'^«llfi* (D) — C'||A„e„||L2 /^-| = CA„. Thus, the partial sums Tfc are members of H|j(Z?); hence, 

if fc < 



L\Tl)-L\Tk) 



I 

E ' 

n=k+l 



en>L|j(_D)-^^(en) 



I 

E 

n=fc+l 



The finiteness of the sum '^n>i^n{'''-i^n)]j (d) thus implies that {L^{Tk)\ is a Cauchy se- 



fc>i 



quence, which by virtue of the completeness of L|,j(Z?) converges to some g* £ 1j\j{D). The 
continuity of L^^ implies that the converge to L^^g* in H^^(_D). As the partial sums converge 
in LIj{D) to t, r = i'^g* e Ht^(D). We have thus proved (1537)) . □ 



We intend to exploit the previous single-domain results in order to say something about the 
multi-domain case. To this end, we define, for i £ [N], the distributional operators Li: {ip £ 
Ll^{D): Vg,^ e [WlliDm ^ P'(D) by 

L.iip) := Mridivq^(A/,Vq^^). 

We also define Li := —Li + I, where / is the identity operator mapping D'{D) onto itself. An 
easily verifiable and important property of these operators is that, as long as their composition 
is well-defined, they commute with respect to their spring index. Hence, we can naturally use 
multi-indices in Nq to denote the repeated application of distinct Li or L^: 



Lfi .- L-^ 



oL 



N ' 



L, 



■.^l{^o...oL^ 



N ' 



where any zero among the Pi is assumed to give rise to the identity operator. For these multi- 
indices we define the function |/3|^ := maxigjjv] A. Now, for derivatives in 2?'(D), the multi-indices 
belong to and come naturally grouped in N length-d sub-multi-indices (one for each factor 
of the Cartesian product D = Di x • • • x Dpf). Thus we define the function Moo i ^ ^o'^ ^ by 

l"loo,l l("l'---'"w)loo,l 



max \ ai 

i<£[N] 



max \ai 

i<£[N] 



that is, the maximal order among the component single-domain multi-indices. 

With this compact notation, we now define the Hilbert spaces (with corresponding norms) 



= 1 



l|2 



'{0) 



E I 

l«U,l<2 



\datp\ 



E II^M^)IIU(D) 



} , (5.40a) 
(5.40b) 



l/3|cx= = l 
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and 

H'm""'^(D) := e L^(D): 3^ e LI,{D), , < 4; 

L^(^) G UliD), = 1; Lpiv) e l2,(D), = 2}, (5.41a) 

ll<^lll*-""'(D) - ll^«<y5||L2 (D) + \\Lpif)\\lliD)+ Y II^/3(¥')IIl^(d) • (5-41b) 

l«U,i<4 l/5U = l l/3U=2 

The following lemma holds. 

Lemma 5.12. Forme {2,4}, H;;;''""'(D) C HT<'"'(D). 



Proof. We recall that, by Lemma |5.3| . ((An, e„))j^gp,N as defined in (IS.lip is a complete set of 
solutions of the M-weighted eigenvalue problem (|5.8I) and that the latter have tensor-product 
structure. Also, by the definitions in (j5.13p and (j5.14L Kj^ ' (D) is the space of L^(D) functions 
whose squared Fourier coefficients, weighted with the coefficients defined by 

N 



T^. J. 

= 1 

have finite sum. 

If (ys £ Hj|^^'™"^(D), one can apply to it each operator Li a total of m/2 times cumulatively and 
the result wiU lie in L^(D); i.e., i(m/2....,m/2) (v') G L^(D). By Parseval's identity. 



00 > 



L 



2 2 ^ m 

(m/2,...,m/2)'/' = Y im/2,. . .,m/2){v) , Sri) ^ ^ J2 11 (-^"^ "^"^ (D) ' 

(5.42) 

where the second equality is justified by the density of C|f (D) functions in (D) , the regularity of 
ip and the Cartesian product form of the domain D and the tensor-product form of the Maxwellian 
weight function M. As the finiteness of the last expression in (|5.42p is exactly the condition for 
membership in H^' ' (D), we obtain the desired result. 

□ 



We recall that Theorem 5^ gives a condition on the parameter of the space Hj^*'"'(D) under 
which it becomes a subspace of the abstract space Bi , which in turn is connected by (14.141) to the 
space Ai of fast convergence of the greedy algorithms (cf . Theorem l4.6l and Theorem l4.7t ). Then, 
from Lemma l5.12l it is apparent that the arguably less abstract space Hj^'""''(D) will be a subspace 
of Bi for a suitable choice of the parameter. We shall now make this statement more precise. 

Theorem 5.13. Let H^* '(D) be defined according to (j5.14p . where d € {2,3}, as elsewhere, is 
the common dimensionality of the Cartesian factors that make up D; then, 

Hr''^(D)cHr(D)c6i. 



Proof. Lemma l5.12l gives that H^™''(D) C '(D). As 4 is greater than 1 + f for both d 
and d — 3, Theorem |5.4| gives that HT'"(D) c Bi. 



□ 



Remark 9. If the hypotheses we have been making throughout this work (i.e., Hvpotheses 1X1 [Bl 
[Cl [D1 and |El are met, nothing in our arguments essentially restricts the results to the physic ally 
relevant case s d = 2 and d — 3. In particular, in the case of d = 1, the combination of Theorem l5.4 
with Lemma l5.12l and the fact that 2 > 1 -1- ^ yields that 

H'M'"'''(D)cHr(D)cSi. 

Sobolev spaces of dominating mixed smoothness akin to H^™"^(D) can also be shown to be sub- 
spaces of the regularity class Bi in the case of the classical Poisson problem studied in [LBLM09] : 
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Find e HJ(D) (with the standard meaning of the Sobolev space Hq(D); i.e., the set of all 
elements of H^(D) that have zero trace on 9D — not a zero-weighted Sobolev space!) such that 

('/',¥'>Hi(D) = (/,¥')l2(d) V(^eHj(D), 

where £) = Di ® ■ ■ ■ ® Dm and each Di, i G [N], is an open interval. The corresponding greedy 
algorithms seek approximations that are linear combinations of (^jgjjy] ^ai-^i) functions. The 
argument of Theorem l5.4l above holds in this case without any change, and so, given that the 
n-th eigenvalue of the corresponding analogue to the partial-domain eigenvalue problem (j5.7p is 
proportional to n^, we have that 



if e l2(d) 



E 




)l2(D) 



< oo 



In this non-degenerate setting it is possible to identify the space on the left-hand side of the above 
expression with 



^(D) n Hi(D) := jy^ e : G L'{D), \a[ 



max a, 

l<i<N 



< 2 



This characterization should be contrasted with the condition for membership in Ai (which is 
identical to Bi in this unweighted setting) derived in ILBLM091 Remark 4], which demands, 
instead, that the true solution belongs to H™(D) n Hj(D), with m > 1 + N/2. In fact the 
characterization given in |LBLM09l Remark 4] can be generalized to the requirement that the 
exact solution belongs to H'"(D) n Hq(D), with m > 1 + Nd/2, when the factor domains are no 
longer one-dimensional but d-dimensional; and, thanks to Theorem l5.5l . such a characterization in 
terms of standard Sobolev spaces (rather than spaces of dominating mixed smoothness) also has 
a counterpart in our degenerate setting. 

An attractive feature of spaces of dominating mixed smoothness is that their regularity index 
is independent of N and such spaces are therefore more naturally suited to (high-dimensional) 
tensor-product settings such as ours^We note in this respect that we conjecture that the reverse 
of the inclusion st ated in Lemma |5.12| also holds, implying equality of the two spaces there — 
just as in Lemma Is.lll for the single-domai n sp aces. We s uspec t that the proof of this would 
involve tensorizing the statements of Lemma 5/7 and Lemma Is. id . However, even if Lemma 5.12 
held with an equality of spaces, there would still be some slack between the discussed mixed 
smoothness levels and the lower bound of the admissible parameter m such that H^'"' (D) C Si. 
The reason is that we have gone about obtaining elliptic regularity results by two integer degrees 
of differentiation at a time. Consequently, we have not defined anything we could label W^j,{Di) 
or H^'™"^(D) with m ^ {2, 4} while being consistent with the definitions we have given for the 
single-spring spaces Il|,j. (D^) in (|5.38p and Hf j ^Di) in (|5.39p . and with the multi-spring spaces 

T2,mix 



H^™(D) in and Hj 



4,miXj^Q-j (|5,4ip . Given the presence of the second-order operators of 
the form Mr-'^div(MiV') and M^^div(MV-) in the definition of these even-indexed spaces, it is 
not immediately clear what a suitable explicit definition of the analogous odd-indexed — let alone 
non-integer-indexed — spaces might be. We shall address this question by using function space 
interpolation. 



We start with the fact that, given two nets of positive weights E^'^ — [a. 



e {1,2}, 



one can show that for 9 £ (0, 1) the (real) {9, 2)-interpolation space between them obeys 



HS'''(D),H^^-'(D) _=H^(D), 



e.2 



where Y. — (dr, 



and (7r 



(^fi^^^ (''^" ^ ^ ^^^^ equivalence of norms (the 

proof is a simple modification of the argument given in [Tar07l Chapter 23]). As, according to the 
definition in (j5.14p . 



and 



r(2+2e) 



-(2) 
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for all 6* e (0, 1) and n in , it follows that 



H 



rp(2e) 

M 



(D) 



Hr(D),Hr(D) 



e,2 



rj,(2 + 2e) 

and 



(D) 



H 



rp(2) 
M 



(D),Hr(D 



with equivalence of norms. Given that the inclusions in Lemma 15.121 are actually continuous 



embeddings and H,^ '(D) = L^(D), it follows that 

V2,mix 



L^M(D),Hr"'(D) 



C and 



with continuous embedding. Since whenever m > 1 + | we have that H,^' '(D) C Si, defining 

H~ m,mix / 
^ 



^(D) as the interpolation space appearing on the left-hand side of the corresponding inclusion 
above if m G (0,2) or to e (2,4) and as L^(D) if m = is an appealing idea, for then we can 
simply state that 

2 



TO > 1 



2 ' ^ 

6. Conclusions and directions for future work 



We proved the well-posedness (Theorem l3.2r ) and convergence fTheorem l4.3l and Theorem l4.5r ) of 
two greedy algorithms, which seek approximations to solutions of high-dimensional and degenerate 
Fokker-Planck equations using a separated representation procedure. We then gave sufficient 
conditions on the true solution of the equation for the fast-convergence of the approxima tion s 
given by those algorithms; first, in terms of sum mability of Fourier coefficients (Theorem \5.4 \. 
and then, in terms of regularity (Theorem IS.lSh . In the process of proving these main results, 
a number of auxiliary results were proved, some of which are of interest in thei r ow n right; e.g., 
function spaces with tensor-product weights inherit compact embedding (Lemma [2^ and density 
(Corollary |4.4| ) properties from the spaces corresponding to the weights that appear as factors of 
the tensor product weight; the existen ce of elliptic regularity results for single-spring degenerate 
elliptic problem s (pa rticularly Lemma IS.lOh : and eigenvalue asymptotics in the same degenerate 
setting (Lemma lc.Sl in Appendix 0)- 

The greedy algorithms described in Section are abstract. They entail obtaining the true 
minima of functionals in nonlinear manifolds embedded in infinite-dimensional function spaces 
(cf. p.3p and p. 41) '). Any practical implementation of the separated representation strategy must 
then introduce a discretization (e.g., by a finite element method or a spectral method) and a 
procedure for the approximation of those minima in the resulting discretized manifolds (e.g., an 
alternating direction scheme, Newton iteration, etc.). The mathematical analysis of the effects 
of the discretization and the use of approximate minimization algorithms on the convergence 
of the greedy algorithms is the subject of ongoing research. On a related note, we are also 
interested in the implementation of the combination of the separated representation strategy 
and the alternating direction scheme described in (I1.14p and (jl.lSp in order to approximate the 
full Fokker-Planck equation (jl.4p . Further up in model complexity is the coupling between the 
full Fokker-Planck equation and the Navier-Stokes equations for the velocity and pressure of an 
incompressible solvent, which is also of interest to us. The Navier-Stokes-Fokker-Planck system 
is a fully coupled macro-micro system, since the configuration probability density function given 
by the Fokker-Planck equation feeds into the Navier-Stokes equations a contribution to the extra- 
stress tensor while the Navier-Stokes velocity field enters in the Fokker-Planck equation (cf. 
[ECAH871 §15.2]). An important property of the full Fokker-Planck equation is that its solution 
is almost everywhere nonnegative and has unit integral over the configuration space D (i.e., it is 
a probability density function) at almost every point in time t and space x if the initial condition 
has those properties. It is of interest to learn whether the separated representation strategy can 
be adapted to give approximations that also preserve the property of being a probability density 
function. 

The generalization of our results to other tensor-product-based high-dimensional PDEs is also 
of interest. In particular, the adaptation of the separated representation strategy to the Fokker- 
Planck equations for the configuration of bead-rod polymer chains (see, e.g., jBCAHSTl §11.3]) is 
of relevance; these models are not covered by our arguments, at least not in their present form. 
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Appendix A. Basic results for Sobolev spaces weighted with CPAIL Maxwellians 

We shall derive some key properties of the function spaces associated with the CPAIL force 
model (jl.3p . with parameter 6 > 3, using the corresponding properties of the function spaces 
associated with the FENE force model (|1.2p . with parameter 26/3. 

Let 6 > 3. It follows from (jl.2p . (jl.3p and (11.51) that the Maxwellian Mq associated to a spring 
obeying the CPAIL model with parameter h and the Maxwellian Mp associated to a spring obeying 
the FENE model with parameter 26/3 are, respectively, 

Mc(p) = Zcexp(-|p|V6) [l-^] , p e = S (o, c K'^ 



and 



Mf(p) = - , p e i^F - i? (o, c 

where Zq and are positive constants whose specific values are of no particular relevance below. 
Let us denote by T the invertible map p G Dq i-> ■\/2/3p G Dy- On defining M : Dq — > R via 
M := Mp o T we find that there exist positive constants ci and C2 such that ciM < Mq < C2M. 
This implies that ll\.j^{Dc) and H^-^(_Dc) (the latter is well-defined since inherits from 

M^^ its ^'loci^c) regularity — thereby falling under the hypotheses of |K084[ Theorem 1.11]) are 
algebraically and topologically the same space. The same is true of the pairs of spaces given by 
Lmc(^c) and Ll^(Dc) and H(Mc;Z?c) and H(M;Dc). 

Now, T and are \C°°(Dc)]'^ and [C°°{Dp)]^ functions, respectively. Then, an argument 
analogous to Lemma [b . 3l leads to the fact that the composition with is a well-defined, invert- 
ible, linear and bounded operator between ]i'-^^{Dc) and 'H\.j_^{Df) and also between L^(I?c) and 
L^p(Df), and its inverse is the composition with T. By (|2.3p . composition with is also such 
an operator between ll{Dc', M) and II(I?f; Mp) having as its inverse the composition with T. 

We can thus use the connection between the A/F-weighted spaces and the M-weighted spaces 
and the connection between the latter and the Afc-weighted spaces to state that 

and 

-H(LIf:A/f) , ^ T 7 -T-H(Dc;Mc) 

C^iDp) "=R{Dr,MF) =^ {/oT: /GCg^pF)} ^'=Hpc;Mc). 

As 2&/3 > 2, the statements on the left-hand side of the above implications hold (as noted in 
Remark!^ and Remark 0); consequently, so do the statements on each right-hand side. By noting 
that, on account of its infinite differentiability, the composition with T maps C^{Df) into Cg°(£'c) 
and that Mc itself is a C°°{Dc) function, we have proved the following lemma. 

Lemma A.l. Let M : D ^ ]S. be the Maxwellian associated to a spring obeying the CPAIL force 
model (|1.3p with parameter 6 > 3. Then, the compact embedding ^{^^{D) ^ h\,j{D) holds; and the 
set G^{D) IS dense in H(D;M). 

Appendix B. Some results on distributions 
Throughout this section 51 will denote an open subset of R'^. 

Lemma B.l. Let a G Nq and let f G Ll^^{n) and g G C(il) be such that dpf G Lj'Qj.(il) and 
dpge C{n) for all (3 < a. Then, dJifg) G hl^^n) and 

dM9)-Y.^-^^0,fO..,g, (B.l) 
where the convention 7! — nje[d] 7 ^ has been followed. 
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Proof. The result is obviously true for a — (0, ... ,0). Then, in the |a| ~ 1 case, the result is 
stated under the assumption of /, daf, dag, fg and fdag + gdaf being members of Lj'Q^(n) 
(which is clearly implied by our hypotheses) in the discussion that follows Theorem 7.4 of [GTOlj . 
The final result follows from standard combinatorial arguments and an induction procedure. 

□ 

The purpose of the following lemma is to formulate a result analogous to Theorem 3.41 of |AF03) 
for weighted Sobolev spaces without resorting to density arguments, which may be unavailable for 
one or both of the weighted Sobolev spaces being connected. 

Lemma B.2. Let T he an invertible C°°(fl) transformation with codomain fl and let f € Lj'^^(ri) 
be such that its distributional derivatives are in Ll^^{n) up to the order a e N'*. Then, 

OMoT-')^ J2 MaAdpfoT-^)^l.l,{h), (B.2) 

1<I/3|<I«I 

where Ma^/s is a polynomial of degree not exceeding \/3\ in derivatives of orders not exceeding \a\ 
of the various components of . 

Proof. Let S denote the inverse of T and let Sk denote the its fc-th component. From Theorem 6.1.2 
in |Hor83l and the remark that follows it we know, first, that there exists a unique continuous 
linear map S* : I?' (17) — > V {Vl) whose restriction to G{Vl) is m H> uo S* and, second, that the chain 
rule, 

d 

d.jS*u^^djSkS*dkU 

k=l 

holds in 2?'(17). It is easy to see (either directly or from the proof of Theorem 6.1.2 of |H6r83j ) 
that S*u has the explicit form 

S*u{ip) ^u{{ipoT)\det{WT)\) Vip e C^in). 

For a regular distribution such as / the above characterization and the change of variable 
formula for integrable functions (see, e.g., |Bog07[ Theorem 3.7.1]) makes S* f precisely the regular 
distribution associated with the Ll^^{fl) function f o S. Similarly, S*dkf will be the regular 
distribution associated with the Lj'Q^(j7) function dkfoS. Hence, dj{foT^^) — J2t=i djSkdkf o S 
and (|B.2p is proved for \a\ = 1. An induction argument then establishes (IB. 21) in the general 
case. □ 



Lemma B.3. Let Q and T be as in Lemma \B.Sl and let w be a weight function defined on fi. 
Then, f £ H™(n) if, and only if, f o T^^ G H™(fJ) and there exist positive constants Ci(m) and 
C2{m) such that 

ciWf °T-^\\-a-^(n) < II/IIh-co) < C2I!/ o r~^||Hn.(o)> 

where w ~ w o T^^ . 



Proof. We use Lemma lB.2l to replace the first part of the proof of Theorem 3.41 of [AF03j . Then, 



the rest of that proof, mutatis mutandis, carries over to our case. □ 

Appendix C. Eigenvalue asymptotics for Ornstein-Uhlenbeck operators with 
fene and cpail potentials via the llouville transformation 

Lemma C.l. Let Vl C be a hounded and convex domain of class and let w G C^(r2) he 
a positive function such that C§(f2) is dense in H^(f2) and Hj^,(f7) (g L^(r2). We further assume 
that 

(1) inf Qi{p) > ~oo, or 

(2) there exists a 9 > such that 70 := inf c)(p)^(3e(p) G (-1/4,0], 

pen 
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where Qq := Q — w ^/^div(wViy ^/^) and is the distance-to-the-boundary function in Q. 

Let (A„)„gpj be the (ordered, with repetitions according to multiplicity) sequence of eigenvalues 
of the problem: Find A G M and u G H^, (fJ) \ {0} such that 

(u,«)hj„(o) = A(u,i;)L.,(i,) V«eH^(r!). (C.l) 

Then, there exist positive numbers ci and C2 and a natural number uq such that 

n>no =^ cin^/'* < A„ < caTi^/'*. (C.2) 

Proof. Let, for > 0, (Ae.n)„gj!j be the (ordered, with repetitions according to multiplicity) 
sequence of of eigenvalues of the shifted problem: Find A® G M and u e Hj„(ri) \ {0} such that 

{u,v)iii(n),0 := (Vu, Vw)[L2 + 0(u,w)lj,(o) = A®(u, i;)l2^(o) Vw e Hi,(f^). (C.3) 

By the hypotheses of t he l emma the existence and the accumulation at oo only of the Ae,n is 
guaranteed via Lemma l5.ll . It further follows from the spectral theory of self-adjoint compact 
operators that Ae,n can be characterized by the Courant-Fischer-Weyl min-max principle: 

Ae n = mm max — ; : = mt sup — ; : , (*-'-4) 

scH,i„(n) scc^(n) 

the second equality being a consequence of the density of Co(il) in Hj„(ri) (cf. |Dav95[ Theorem 
4.5.3]). Note that when = 1 the problem (jC.3p and the problem (|C.1I) coincide (and so do the 
sequences (Ae,„)„gN and (A„)„gpj). 

Let L := w^^^^ € C^(0), let z be an arbitrary Cq{Q) function and let y := 



j^{\V{Lyf + Q{Lvf)L- 

[ \Vy\'+ f (e + L-'\WL\')y^+ [ L-'VL-V{y^) 

I |Vy|'+ / (e + i-2|VL|2)2/2_ [ div {L-'VL)y^ 
Jn Jn ^ ' Jn 

[ \yy\'+ [ [e-LdW{L-'\/L)]y' 
Jn Jn 



n 



Similarly, |l2:||^2 (^^-^ = l|y|lL2(o)- As z e Cq{D,) is arbitrary and z i-^ L -'^z is a bijection of CQ(r2) 
into itself, (jC.4p begets 

_ \\^y\\'vHn)]-+InQBy^ 

^e,n — .. ini sup 



If condition ([l} holds, there must exist a 8 > such that Qq > in £7. For such a 9, of course, 
/n Q&y^ ^ 0- On the other hand, if condition ^ is met, then with the particular G given in the 
condition we have that 



the last inequality being a multi-dimensional Hardy inequality (see, e.g., [MMP981 Theorem 11], 
bearing in mind that 7e has been assumed to be nonpositive) . In either case, we can write 

Ae.n> mf sup "„ J" , (C.5) 

where 

1 if condition ([T]) holds. 



< a := 



(1 + 7e/4) if condition © holds. 
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The C'^ regularity of dfl implies the existence of an Eq E (0, 1) such that for each e E (0,eo) there 
exists a subdomain ft^ ^l that is also of class and has measure (1 — e) Fixing e E (0, Eq), 
the fact that the extensions by zero of functions in Co(ile) form a subspace of CQ(ri) and (IC.4P 
imply that the eigenvalues of the unshifted problem (jC.ll) can be bounded from above according 
to 

A„ < mf sup r— . 

dim(S)=ri 2g5\{o} (^;^/L^(n,) 

Now, the right-hand side of (jC.5l) and the right-hand side of (jC.6p are precisely the n-th 
eigenvalue associated with the (variational form of the) problem 

—aAy = iiy in n, y = on dfl 

and the problem 

— div(wVy) + wy — vwy in VL^, y — Q on dfl^, 

respectively. These standard eigenvalue problems obey Weyl's law (this results from the fairly 
general Theorem 2.4 of |Cla67] with input from the regularity result in |Bro61[ Theorem 2.4] — 
alternatively, see |CH531 §VI.4.4]); that is, 

^^oo jyrf/a (2V7F)''r(l + d/2) V / : \ i 

where C := ((2-y/7r)'^r(l + . Particularizing these limits to fi = fin and h' — i^n they 

turn into statements about the rate of growth of the eigenvalues themselves, as opposed to the 
counting functions. That is, 

lim Mn/"'/'' = and lim i^^/n^/^ ^ {1 - e)-^I^C~^l^. 

n— >oo n— ^C30 

From the definition of the shifted eigenvalue problem (jC.3[) . for any 0, it is immediate that 
Ae,ri = A„ + — 1 for all n g N. We then deduce, via the inequalities (|C.5I) and (|C.6p . that the 
asymptotic bounds (|C.2I) hold. □ 

Remark 10. 



(1) It follows from the proof of Lemma lC.ll that, if condition ([!]) holds, the constants ci and 
C2 of (jC.2l) can be taken arbitrarily close to C^^'^ and, consequently, to e ach other. 

(2) One might relax the condition of convexity of the domain in Lemma IC.ll at the possible 
cost of having a stricter lower bound for 79 in condition ^ , as the constant for the Hardy 
inequality might deteriorate. The regularity condition on the domain can be drastically 
relaxed (see, for example |BS70p : however, the literature tends to force one to choose at 
most two among readability, the size of the class of problems covered, and frugality in 



terms of hypotheses. For our purposes, the statement in Lemma lC.ll suffices 



Corollary C.2. The eigenvalues of the eigenvalue problem (j5.7D associated with both the FENE 
model (jl.2p and the CPAIL model (|1.3p obey \C2\ if their parameter hi is greater than 2 and 3, 
respectively. 



Proof. We shall apply Lemma IC.ll . For both the FENE and CPAIL models the domains (being 
balls) and their associated Maxwellian weights are regular enough. The compact embedding and 
density hypotheses are satisfied in the parameter ranges under consideration (cf. Hypothesis [b|. 
Remark y. Remark iand It only remains to prove condition ([T]) or condition ([2]). 

From (|1.2p and (|1.5p it follows that the Maxwellian associated to the FENE potential is 



M,{p)^Zr^ {l~\p\yb,f^\ PEB{0,^), (C.8) 
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where Zi is a positive constant. A direct calculation returns that with this weight Qe is 

«e(.,^e.(i-i)„r(i-M!)-\ 

In this form, it is readily apparent that Q\ is bounded from below in its domain -8(0, \fbi) (i.e., 
(H]) holds) if hi > 4. From the fact that d{p) = Vb- \p\ for all p in the domain under consideration 
it is easy to see that d^Qs is always bounded from below and uniformly continuous up to the 
boundary. If bi € (2, 4], Qq is never bounded from below, so it takes negative values and thus the 
infimum of d^Qs is strictly less than zero. As 0^ is continuous and positive within the domain 
yet zero at its boundary, the existence of a that makes case ([2]) hold is equivalent to demanding 
that 

lim 5(p)2gi(p)e (-1/4,0]. 

As in the range hi £ (2, 4] that limit is bi{bi/A — l)/4 we see that the condition ^ holds there. 
Analogously, (jl.3D and (11.51) imply that the Maxwellian associated to the CPAIL potential is 

M,(p) = Z-iexp(-lpl76) (1-Ipl7b.)''/', peB(0,v/^), (C.9) 
with Zi a positive constant. Again, a direct calculation yields 

By arguments similar to those given when considering the FENE potential, we have that condition 
([T]) holds if 6j > 6 or if 6,; = 6 and d = 2; and that condition ^ holds if 6i G (3, 6]. □ 

If two weights w and w defined on a domain are comparable — that is, there exist two positive 
constants ci and C2 such that ciw < w < C2W — a number of consequences follow immediately. As 
discussed elsewhere, L^(r2) and L'^{il,) on the one hand and Hj„(ri) and Il)j(r2) on the other will 




be one and the same algebraically and topologically. In particular, the hypotheses of Lemma 1 5.1 
will be met by the eigenvalue problem 

(e,w)Hj„(0) = A(e,'u)L2 (fj) Vu G Hi,(fi) 

if, and only if, they are met by the eigenvalue problem 

(e,w)H^(n) = ^(e,w)L|(o) Vu € H^^j^^. 

The inf-sup characterization (cf. (jC.4p ') of the successive eigenvalues of both problems allow for 
the bounds 

— A„ < A„ < — A„. 

C2 ci 

That is, the bounds (|C.2p will hold for one set of eigenvalues if, and only if, they hold for the 
other. This allows for establishing the following sufficiency condition for weights defined on two- 
or three-dimensional balls, which is in most cases much easier to test than the conditions of 



Lemma C.l 



Lemma C.3. Let be an open ball in two or three dimensions and let w be a positive and 
continuous weight defined on f2 with the property 

for all p £ ^ such that c)(p) < 5, for some exponent a > 1, for some margin 5 > and some 
positive constants ui and a-i . 

Then, the eigenvalues of the problem 

(e,w)Hi,(f2) = A(e,?;)L2^(o) Vw G Hi,(f^) 

obey the two-sided bounds (jC.2p . 
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Proof. If the radius of the bah happens to be v 2a the conditions on w force it to be comparable to 
the FENE MaxweUian (|C.8p and so the result follows from the above discussion. Otherwise, one 
just needs to rescale the domain; this will effect a fixed linear transformation on the eigenvalues, 
but will not affect the validity of the bounds (|C.2|) (the constants involved will change, though). □ 



Remark 11. The eigenvalue problem (j5.7l) associated with either the FENE or the CPAIL model 
falls within what is called weak degeneracy case in the Russian spectral theory literature; i.e., 
problems of the form: Given fl C R'*, find (A, u) e R x (H^c (n) \ {0}) such that 

/ (AS/u-S/v + huv)^"' ^ \ I buvd/^ Vi;eHj„(f7), (C.IO) 
Jn Jn 

where a — /? < 2/d (see [VS741 §1] for the precise statement, which includes additional conditions 

on n, A, h, b, a and /]). As, in the FENE and CPAIL versions of (|5.7I) . the same weight (the 

associated MaxweUian) appears in both the left- and right-hand side bilinear forms, and, in both 

cases, that weight is bounded from above and below by powers of c) (cf. (jC.Sp . (|C.9|) '). it turns out 

that our problem is equivalent to a problem of the form (|C.10p with a — f3 = 0. 

The result, according to |VS74[ Theorem 1.1] and assuming that 6 > is that 

lim A-'^/2^{n e N: A„ < A} = , , — / , (C.ll) 

^ ^ (2V^)''r(l + d/2) v/dit(l) ^ ' 

(compare this with (|C.7I) : note also that in |VS74] the statement is made in terms of what in 
our notation is 1/A). The problem with this particular source is that, for a proof, it remits the 
reader to either one of two publications. The first, |BS72) proves related yet not directly applicable 
results — there is a gap that needs to be bridged by means, perhaps elementary, that are unknown 
to us. We have not been able to get hold of the second, |Tas75] by G. M. Tascijan (also romanized 
as Tashchiyan). However, the latter is also cited in [TasSll Theorem 1], where a generalization of 
(IC.11|) is proved, under the condition (in our notation) d > 2. 

Acknowledgement. We are grateful to Professor Marco Marietta (Cardiff University) for 
helpful suggestions regarding the Liouville transformation. 
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